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1.  INTRODUCTION 

1.1.  Background 

Motor  vehicle  pollution  comes  from  the  by-products  of  the  imperfect  fuel 
combustion  process  (exhaust)  in  the  engine  and  from  the  evaporation  of  the  fuel  itself. 
The  main  pollutants  are  hydrocarbons  (HCs),  nitrogen  oxides  (NOJ,  and  carbon 
monoxide  (CO).  HCs  exhaust  emissions  result  when  fuel  molecules  in  the  engine  do  not 
burn  or  burn  partially.  HCs  also  escape  into  the  air  through  fuel  evaporation  that  occurs 
in  several  ways.  The  diurnal  evaporation  resulting  from  the  fuel  tank  heating  as  the 
temperature  rises  during  the  day,  the  running  losses  of  the  hot  engine  and  exhaust  system, 
the  hot  soak  after  the  vehicle  is  turned  off,  and  the  refueling  losses.  NOx  exhaust 
emissions  are  formed  under  high  pressure  and  temperature  in  the  engine  when  nitrogen 
(N)  and  oxygen  (O)  atoms  in  the  air  react  to  form  various  NOx.  CO  is  a  product  of  the 
incomplete  combustion  that  occurs  when  carbon  in  the  fuel  is  partially  oxidized  rather 
than  fully  oxidized  to  carbon  dioxide  (C02).  The  chemistry  of  the  formation  of  these 
chemicals  is  simplistically  represented  as  follows:  Fuel  (HCs)  +  Air  (02,  NJ  -> 
Unburned  HCs  +  CO  +  NOx  +  C02  +  H20.  02,  N2,  and  H20  represent  the  molecules  of 
oxygen,  nitrogen,  and  water  respectively.  In  this  research  project,  the  exhaust  CO  and 
HCs  emissions  are  the  main  focus.  Substantial  amounts  of  CO  are  produced  whenever  a 
rich  mixture  is  used.  Spark-ignition  engines  employ  rich  mixtures  during  startup  to 
prevent  stalling  and  at  wide-open  throttle  conditions  (hard  acceleration)  to  provide 
maximum  power.  For  stochiometric  and  lean  mixtures,  CO  is  produced  in  large  quantity 
at  typical  combustion  temperature  as  a  result  of  the  dissociation  of  C02.  Partial  oxidation 
of  oil,  and  deposit  layers  of  unburned  fuel  is  another  major  source  of  CO,  especially  in 
lean  operation.  Unburned  HCs  may  result  from  the  flame  quenching  process  which 
leaves  thin  layers  of  unburned  fuel-air  mixture  at  the  cylinder  walls  and  at  the  entrance  to 
crevices  such  as  piston  top  lands  and  ring  packs.  Adsorption  and  desorption  of  fuel  into 
oil  layers  in  the  cylinder  walls  are  also  great  contributors  to  unburned  HCs.  For  lean 
and/or  dilute  mixtures  approaching  the  flammability  limits,  unburned  HCs  may  result 
from  incomplete  flame  propagation  in  the  bulk  of  the  charge  (Turns,  1996).  The  effects 
of  CO  and  HCs  on  the  environment  and  human  health  are  well  documented.  They  alter 


the  properties  of  the  atmosphere  and  precipitation,  harm  vegetation,  deteriorate  materials, 
and  increase  morbidity  and  mortality  in  humans  (Seinfeld,  1986).  For  example,  HCs 
react  in  the  presence  of  NOx  and  sunlight  to  form  ground-level  ozone  (03).  03  irritates 
the  eyes,  damages  the  lungs,  and  aggravates  respiratory  problems.  Moreover,  many 
exhaust  hydrocarbons  are  carcinogens.  CO  is  a  colorless  and  odorless  chemical  known  to 
reduce  the  flow  of  oxygen  in  the  bloodstream. 

Construction  work  zones  greatly  affect  highway  traffic  patterns  and  vehicular 
exhaust  emissions.  As  a  matter  of  fact,  the  reduction  in  the  number  of  driving  lanes 
causes  heavy  congestions  and  long  delays.  Because  of  the  unstable  stop-and-go  effect  of 
traffic  flow,  the  high  density  of  vehicles,  and  the  variable  fleet  mix,  very  large  amount  of 
exhaust  emissions  are  generated;  sometimes  exceeding  the  national  standard  levels.  The 
current  federal  standards  for  1994  models  and  later  expressed  in  gram  per  mile  (g/mile) 
for  exhaust  emission  from  cars  are  0.25  HCs,  0.4  NOx,  and  3.4  CO.  Despite  the 
implementation  of  the  standards  and  the  development  of  efficient  exhaust  emission 
control  technologies,  the  improvement  of  air  quality  has  not  been  great.  This  is  due  to 
three  important  factors.  First,  the  slow  introduction  of  new,  cleaner  cars  due  to  the 
increase  of  the  average  age  of  the  vehicle  fleet.  Second,  the  offset  of  the  emissions  gains 
by  the  growth  of  vehicle  miles  travelled  (VMT)  which  is  doubling  every  twenty  years. 
Third,  and  perhaps  the  most  important  is  the  inadequate  understanding  of  real-world 
emissions.  This  results  in  emission  test  procedures  and  standards  that  do  not  represent 
the  actual  source  of  emissions. 

Current  mobile  source  emission  estimation  is  the  product  of  two  steps.  The 
determination  of  a  set  of  emission  factors,  and  the  estimation  of  vehicle  activity.  The 
emission  factors  are  produced  based  on  average  driving  characteristics  of  the  federal  test 
procedure  (FTP)  which  has  a  specified  driving  trace  of  speed  versus  time  intended  to 
reflect  actual  driving  conditions  on  both  arterial  and  highway.  Emissions  of  CO,  HCs, 
and  NOx,  are  collected  for  three  cycles  (cold  start,  hot  stabilized,  and  hot  soak)  and  are 
used  as  base  emission  rates.  Then  correction  factors  are  applied  to  adjust  the  basic 
emission  rates  to  reflect  the  observed  differences  between  different  modes  of  operation. 
Temperature  and  speed  correction  factors  are  also  used  to  adjust  the  emission  rates  of 
non-FTP  cycles;  speed  greater  than  57  mph,  and  acceleration  greater  than  3.3  mph/s. 


Vehicle  activity  including  VMT,  speed  distribution,  number  of  vehicles,  number  of  trips, 
and  vehicle  mix  are  derived  from  macroscale  transportation  models.  This  emission 
inventory  method  has  several  shortcomings  including  the  inaccurate  representation  of 
actual  driving  behavior,  and  the  flaw  in  the  emission  factor  estimation.  Because  of  the 
characteristics  of  this  methodology,  it  cannot  be  used  efficiently  to  evaluate  traffic  flow 
improvements  that  are  microscale  in  nature. 

A  Fourier  transform  infrared  (FTIR)  air  monitoring  system  can  measure  real-time 
on-road  exhaust  emission  concentrations  from  vehicles  operating  in  various  driving 
modes  instantaneously.  Such  a  system  consists  of  an  infrared  source  which  transmits  a 
collimated  beam  that  is  detected  by  the  spectrometer  which  contains  the  interfering  optics 
and  the  detector.  Post-car  signals  immediately  acquired  after  a  vehicle  passes  through  the 
infrared  beam  are  compared  to  the  background  ambient  air  signal  in  order  to  determine 
the  concentration  of  the  pollutant  in  the  exhaust  plume.  Since  the  vehicle  fuel  economy 
cannot  be  determined  in  real-time,  the  concentration  measured  by  FTER.  systems  are 
expressed  on  a  volumetric  basis  (i.e.  ppm,  %).  The  rapid  introduction  of  this  monitoring 
system  has  prompted  some  states  to  develop  new  set  of  standards  based  on  FTIR  data. 
For  example,  The  North  Carolina  current  standards  for  198 land  later  models  of  light  and 
medium  duty  vehicles  are  12000  ppm  CO  and  220  ppm  HCs.  The  California  HCs 
standard  for  1994  models  of  light  and  medium  duty  vehicles  is  200  ppm  HCs. 

1.2.  Objectives 

The  main  objective  of  this  research  is  to  develop  a  real-time  modal  exhaust 
emission  concentration  model  for  light  and  medium  duty  vehicles  using  real-time 
emission  concentration,  traffic,  and  meteorological  data.  The  data  were  collected  during 
the  Purdue  Vehicle  Emission  Monitoring  and  Modeling  Project  (PVEMP)  field 
experiments  conducted  in  the  spring  and  summer  1997  on  the  Borman  Expressway.  The 
Borman  Expressway  is  a  multi-lane,  16-mile  link  of  highway  1-80/94  located  in  northwest 
Indiana  (Figure  1.1),  a  non-attainment  area  for  03  and  other  criteria  pollutants.  The 
traffic  volume  on  the  Borman  Expressway  is  extremely  high  and  alternative  routes  are 
limited.  Consequently,  transportation-related  emissions  cause  major  air  quality  problems. 
The  implementation  of  the  model  and  the  results  of  the  project  will  integrate  air  quality  as 


another  parameter  in  the  traffic  flow  improvement  of  the  Borman  Expressway,  especially 
in  construction  zones.  The  operation  of  the  model  will  require  the  use  of  real-time  traffic 
parameter  inputs  obtainable  only  from  advanced  traffic  monitoring  systems  such  as  the 
Autoscope™  or  similar  systems.  This  paper  presents  the  results  of  the  monitoring 
program,  the  data  reduction  and  analysis  effort,  the  development  of  the  real-time  modal 
exhaust  emission  concentration  model,  and  the  implementation  suggestions  of  the 
research. 
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Figure  1.1.  Local  area  map  showing  the  Borman  Expressway  (source:  INDOT  web  page). 


2.  FTIR  REMOTE  SENSING  TECHNOLOGY 

The  FTIR  technology  is  used  for  the  detection  of  chemical  compounds  through  a 
particular  medium.  An  energy  source  in  the  form  of  infrared  (IR)  radiation  is  sent 
through  a  desired  medium.  This  radiation  causes  the  atoms  in  a  molecule  to  resonate  at  a 
certain  frequency  specific  to  the  molecules  mass  and  structure.  When  the  particular 
frequency  in  the  IR  band  meets  a  molecule  resonating  at  the  same  frequency,  the  IR 
radiation  is  absorbed.  This  causes  a  reduced  energy  level  at  that  specific  wavelength  to 
be  measured  by  the  sensing  device.  Previous  tests  have  been  conducted  and  reference 
libraries  have  been  established  to  aid  in  determining  the  wavelength  at  which  the 
molecules  of  concern  absorb.  These  libraries  are  referenced  by  the  computer  and 
matched  with  the  output  from  the  spectrometer  to  ultimately  identify  the  compound 
existing  in  the  media  being  tested. 

Often  the  FTIR  spectrometer  is  found  in  a  single  unit  on  a  laboratory  bench  and 
includes  a  closed  cell  containing  the  sample  to  be  analyzed.  In  order  to  sample  vehicle 
exhaust  emission  concentrations  in  real  time  however,  an  alternative  setup  has  been 
designed.  Such  a  setup  is  referred  to  as  an  open-path  FTIR  sensing  device  and  is 
identified  by  the  separate  IR  source  and  spectrometer  modules.  The  free  space  between 
the  source  and  the  receiving  optics  is  confined  only  by  the  cylindrical  wave  pattern  of  the 
IR  source.  This  system  includes  the  IR  source,  the  transmitting  optics,  the  receiving 
optics,  the  interferometer,  and  the  detector.  A  basic  macroscopic  schematic  is  illustrated 
in  Figure  2.1. 

The  system  used  for  this  study  is  designed  by  Midac  Corporation.  The  IR  source 
is  a  silica  carbide  (SiC)  element  that  glows  at  a  temperature  of  1550  °K.  The  element  is 
mounted  at  the  center  of  a  twenty-inch  parabolic  reflector  that  yields  a  cylindrical  IR 
wave  pattern  for  easier  alignment  with  the  spectrometer.  Within  the  spectrometer  are 
housed  an  interferometer,  a  beamsplitter,  and  a  detector.  A  Plane  Mirror  Michelson 
Interferometer  is  used  in  the  Midac  spectrometer  (Figure  2.2).  It  consists  of  two 
perpendicular  mirrors  and  a  zinc  selenium  (ZnSe)  beamsplitter. 
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Figure  2.1.  Basic  schematic  of  an  open-path  FTIR  system  (Source  EPA,  92). 
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spectrometer  (Source:  Midac  Web  page). 


After  the  ER  radiation  passes  through  the  medium  to  be  sampled,  it  enters  the 
spectrometer  where  the  beam  is  split.  Some  of  the  radiation  is  reflected  to  a  stationary 
mirror  then  reflected  directly  into  the  detector.  The  waves  that  pass  through  the 
beamsplitter  will  reflect  off  a  moving,  laser-controlled  mirror.  By  moving  back  and  forth 
at  a  high  frequency,  this  mirror  reflects  the  full  range  of  wavelengths  to  the  detector.  A 
basic  schematic  of  this  system  can  be  seen  in  Figure  2.3. 


Moving 
Mirror 


Detector 


Figure  2.3.  Components  of  a  Plane  Mirror  Michelson  spectrometer  (Source:  EPA,  92). 

The  two  mirrors  serve  to  align  the  two  separate  IR  waves  such  that  they  meet  at  the 
detector  with  some  combined  amplitude.  The  energy  absorbed  by  the  medium  is  detected 
by  the  change  in  amplitude  of  the  IR  wave  at  the  detector.  This  concept  is  illustrated  in 
Figure  2.4.  The  detector  is  a  liquid  nitrogen  cooled  mercury,  cadmium,  and  telluride 
(MCT)  detector  that  generates  an  interferogram  in  the  form  of  voltage  at  the  path 
difference  based  on  the  location  of  the  two  mirrors. 
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Figure  2.4.  Path  difference  fundamental  diagram  (Source:  EPA,  92). 


A  computer  and  software  are  used  to  further  analyze  the  interferogram  and 
ultimately  generates  a  concentration  reading  for  the  compound  of  interest.  The  Midac 
system  comes  with  a  spectral  analysis  software  package  called  AutoQuant™;  a  user- 
friendly  version  of  Grams32™.  The  interferogram  generated  by  the  detector  is  reduced 
with  a  fast  fourier  transform  to  generate  a  single-beam  spectrum.  A  single-beam 
spectrum  is  a  plot  of  inverse  interferogram  intensity  versus  the  wavenumber  (inverse 


wavelength).  This  spectrum  is  compared  to  a  previously  recorded  background  spectrum 
to  generate  a  transmission  spectrum.  An  absorbance  spectrum  is  then  calculated  and  is 
defined  as  the  logarithm  of  the  transmission  spectrum.  This  absorbance  spectrum  is  a 
data  set  of  absorbance  versus  wavenumber.  A  schematic  of  the  process  to  this  point  is 
illustrated  in  Figure  2.5.  A  predetermined  absorption  band  is  set  for  each  compound  of 
interest  and  is  further  analyzed  to  determine  a  concentration.  Beer's  law  is  used  and  is  as 
follows:  A  =  (a)CA.  where,  A  is  the  absorbance,  a  the  absorption  coefficient,  C  the 
concentration,  and  A  the  path  length.  The  absorbance  data  is  referenced  to  a  data  library 
of  known  concentration,  path  length,  and  absorbance  spectrum.  The  overall  equation  is 
as  follows: 


^unk  ^  unk  ^  unk 

Where  the  subscripts  'ref  and  "unk"  refer  to  the  reference  values  and  unknown  values 
respectively.  Solving  this  equation  for  C^  generates  the  equation  as  follows: 

unk  ~       2       A 

^unk  ^ref 

This  is  the  last  step  in  a  single-component  analysis.  For  this  study,  both  HCs  and  CO  are 
studied.  For  this  multicomponent  analysis,  a  classical  least  squares  (CLS)  algorithm  uses 
a  ratio  of  path  lengths  to  define  the  multiplier  used  to  determine  the  unknown 
concentration.  A  classical  least  squares  fitting  algorithm  simply  performs  a  linear 
regression  over  a  wave  number  region,  using  the  unknown  and  reference  spectra. 
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Figure  2.5.  The  software  function  diagram  (source:  EPA  92). 
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3.  FIELD  EXPERIMENTS 


3.1.  Preliminary  Study 

The  initial  task  was  to  install  the  16-bit  interface  direct  access  memory  (DMA) 
card  into  the  486DX2-66  computer  that  would  be  used  for  the  research.  The  computer 
was  configured  accordingly  and  the  next  step  was  to  confirm  that  the  spectrometer  was 
communicating  with  the  computer.  To  accomplish  this,  the  system  was  set  up  in  the 
laboratory.  The  appropriate  reference  files  were  located,  installed,  and  referenced  within 
the  software.  Once  a  background  scan  was  taken,  a  bottle  of  hexane  was  opened  in  the 
ER  path  to  confirm  detector  response  (Figure  3.1). 


Figure  3.1.  Preliminary  laboratory  setup  with  hexane  bottle  in  place. 


After  the  laboratory  test  was  complete,  a  series  of  trial  runs  were  conducted  in  a  parking 
lot,  across  a  local  road,  and  in  a  construction  zone  on  1-65  south  of  Lafayette.  The  first 
trial  runs  were  conducted  with  stationary  emissions  sources  (Figure  3.2). 
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Figure  3.2.  A  stationary  trial  run  using  a  1991  Suzuki  Sidekick. 


At  this  stage  in  the  project,  time  was  taken  to  become  proficient  with  the  alignment  and 
calibration  of  the  instrument  and  also  to  identify  the  height  at  which  tailpipe  emissions 
were  best  detected.  Once  calibrated  and  aligned,  a  background  scan  was  taken  and  then 
the  vehicles  were  started.  After  the  vehicles  had  run  at  idle  for  a  period  of  time,  the 
engine  was  revved  and  the  results  observed.  Figure  3.3  shows  the  HCs  exhaust  emission 
concentration  time  series  of  the  stationary  trial  run  performed  on  four  vehicles  on  March 
20,  1997  at  the  parking  lot  behind  the  Purdue  University  Civil  Engineering  building. 
There  are  two  alignment  procedures  that  must  be  followed  when  using  the  system.  First, 
the  alignment  of  the  IR  source  with  the  target  angle  of  the  spectrometer.  This  is 
accomplished  using  a  scope  on  the  IR  source  and  adjusting  the  height  and  angle  so  that 
the  source  reflects  in  the  receiving  optics  of  the  spectrometer  unit.  It  is  also  important  to 
adjust  the  height  of  the  unit  so  that  the  path  that  is  being  detected  is  adjacent  to  the 
tailpipes  on  the  passing  cars.  The  second  alignment  procedure  has  to  do  with  aligning  the 
mirrors  within  the  interferometer  to  record  the  most  efficient  and  accurate  detection.  This 
is  done  with  two  sets  of  screws  on  the  side  of  the  spectrometer  unit.  This  is  a  critical 
alignment  procedure  where  the  parameters  indicated  on  the  computer  are  monitored  as 
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the  mirrors  are  being  adjusted.  Any  vibration  or  bouncing  caused  by  the  passing  cars  has 
an  obvious  effect  on  this  alignment. 
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Figure  3.3.  HCs  Exhaust  emission  concentrations  recorded  during  a  stationary  trial  run. 


A  series  of  dynamic  trial  runs  followed  the  stationary  runs  described  above.  The 
first  was  a  simple  set  up  where  an  occasional  car  would  pass.  Adjustments  in  procedure 
were  made  to  achieve  the  best  detection.  The  second  was  done  on  a  local  four-lane 
roadway  with  a  considerable  amount  of  traffic  (Figure  3.4).  A  third  one  was  performed 
in  a  construction  zone  on  1-65  south  of  Lafayette.  A  considerable  amount  of  data  was 
collected  during  these  trial  runs  and  a  camera  was  incorporated  to  help  identify  the 
vehicle  that  correlated  with  the  emissions  reading.  Figure  3.5  shows  the  CO  exhaust 
emission  concentration  time  series  measured  during  the  dynamic  trial  run  performed  at 
the  1-65  site  on  April  18,  1997.  As  it  can  be  seen,  the  CO  concentrations  were  very 
variable  due  to  the  differences  in  vehicle  types,  age,  and  operating  modes. 
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Figure  3.4:  rour-lane  dynamic  trial  run  with  a  video  camera  set  up. 
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Figure  3.5.  CO  Exhaust  emission  concentrations  recorded  during  a  dynamic  trial  run. 


An  aerodynamic  study  of  the  turbulent  flow  around  and  behind  a  moving  vehicle 
was  conducted  to  help  identify  the  height  from  the  road  that  yielded  the  best  exhaust 
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emission  concentration  measurements  in  the  mixed  volume  behind  a  moving  vehicle. 
Three  basic  vehicle  types  were  identified.  They  included  the  sedan,  the  fastback,  and  the 
boxier  van  type.  Smoke  was  used  to  observe  the  mixing  region  behind  the  van  type 
vehicle  (Figure  3.6)  while  aerodynamic  theory  was  used  to  define  the  other  two  vehicle 
types.  It  was  found  that  a  vigorous  turbulent  cavity  flow  was  attached  to  the  back  of  the 
vehicles  as  they  move.  This  mixed  volume  was  identified  as  the  target  for  exhaust 
emission  concentration  measurements. 
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Figure  3.6.  Aerodynamic  study  of  the  mixing  volume  behind  a  moving  van. 

3.2.  Borman  Expressway  Field  Experiments 

For  the  actual  data  collection  stage  of  this  project,  site  selection  was  of  utmost 
importance  in  order  to  accomplish  the  desired  results.  In  order  to  identify  single  vehicles, 
it  was  necessary  to  limit  the  sampling  site  to  a  single  lane.  Since  the  Borman  Expressway 
is  a  major  three  to  four  lane  corridor  with  very  high  volume  traffic,  it  was  unreasonable  to 
reduce  the  available  lanes  to  a  single  lane.  This  would  produce  not  only  an  incredible 
backup  of  traffic  but  dangerous  operating  conditions  for  the  researchers  as  well.  The 
second  criteria  for  site  selection  was  to  have  a  region  of  likely  acceleration  and 
deceleration  similar  to  what  would  be  characteristic  of  the  stop  and  go  traffic  of  a 
construction  zone.  The  location  that  best  suited  our  needs  was  the  end  of  the  on-ramp 
from  1-65  northbound  to  the  Borman  Expressway  (Figure  3.7). 
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Figure  3.7.  Sampling  site 


ig 


on-ramp  to  the  Borman  Expressway. 


Two  lanes  take  a  reduced  speed,  90°  turn  where  they  merge  into  a  single  lane  before 
further  merging  with  the  55+mph  traffic  of  the  Borman  Expressway.  At  this  location  is 
an  overpass  bridge  just  prior  to  the  merging  of  the  on-ramp  traffic  with  the  main  flow  of 
traffic  on  the  Borman  Expressway.  This  bridge  served  additional  purposes  of  hiding  the 
equipment  from  the  drivers  and  as  a  protective  measure  for  the  equipment  and  personnel. 
The  equipment  was  set  up  just  downstream  of  the  bridge  abutments  on  either  side  of  the 
ramp  (Figure  3.8).  Because  the  equipment  could  not  be  seen  by  drivers  until  it  was 
almost  adjacent  to  the  vehicles,  the  psychological  influence  of  the  equipment  did  not  play 
a  role  in  altering  the  driving  behavior  of  commuters. 

Toward  the  end  of  the  study,  a  location  on  1-65,  two  miles  north  of  the  Frankfort 
interchange  in  the  northbound  lane,  was  selected  to  conclude  the  data  collection  portion 
of  the  study.  This  site  was  downstream  of  a  single  lane  construction  zone  providing  the 
acceleration  that  was  required  (Figure  3.9). 
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Figure  3.8.  Equipment  setup  at  the  Borman  Expressway  sampling  site. 


igure  3.9.  Monitoring  at  1-65  sampling  site. 
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Additional  parameters  were  recorded  simultaneously  with  the  emissions  data. 
These  parameters  included  vehicle  speed  and  acceleration,  video  imageries  for  vehicle 
identification,  and  meteorology  (wind  speed,  wind  direction,  ambient  temperature,  dew 
point  temperature,  and  relative  humidity).  A  camera  was  set  up  under  the  bridge  to 
record  the  images  of  vehicles  as  they  passed  by  the  spectrometer.  Notice  the  tripod  leg 
that  can  be  seen  in  the  dark  left  side  of  Figure  3.8.  This  camera  had  a  digital  time  stamp 
on  the  film  that  was  calibrated  with  the  clock  on  the  computer.  Later,  this  video  data  was 
used  to  generate  a  database  of  vehicle  type  versus  time  at  the  spectrometer.  Speed  was 
measured  with  a  hand  held  laser  speed  gun  located  about  300  feet  upstream  of  the 
spectrometer.  The  speed  data  was  collected  twice  for  each  vehicle  so  that  acceleration 
could  be  determined.  A  small  Hewlet  Packard  (HP)  data  logger  that  through  a  cable 
automatically  recorded  the  speed  data  in  a  form  that  could  be  directly  downloaded  and 
analyzed  handled  the  laser  gun  data  collection.  Like  the  video  data,  the  time  for  the 
speed  gun  was  correlated  with  the  computer  time  for  data  reduction  purposes.  The  third 
data  point  recorded  by  the  laser  gun  was  the  range  to  target  which  would  be  used  to 
determine  the  time  at  which  each  vehicle  passed  the  spectrometer.  This  was  important 
because  congestion,  speed,  and  laser  reading  duration  made  it  extremely  difficult  to  catch 
each  vehicle  when  they  passed  the  spectrometer.  Meteorological  data  were  obtained 
from  the  Gary  and  Purdue  University  airport  weather  databases. 
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4.  DATA  REDUCTION  AND  ANALYSIS 


4.1.  Data  Reduction 

During  the  exhaust  emission  concentration  and  traffic  monitoring,  16,870  vehicles 
were  videotaped.  The  tapes  were  manually  and  visually  analyzed  in  order  to  determine 
the  vehicle  types  as  well  as  correlate  vehicle  occurrence  to  the  time  of  emission 
concentration  measurements.  In  doing  so,  the  correlation  times  included  a  four-second 
time  lag  representing  the  maximum  time  for  vehicle  exhaust  plume  dissipation  (Caddie  & 
Stephens,  1994).  Three  vehicle  types  were  determined.  All  the  automobiles  were 
represented  by  type  I  vehicles.  Type  II  represented  all  the  medium-duty  vehicles 
including  light  and  medium-duty  trucks  and  vans.  Heavy-duty  trucks  were  represented 
by  type  III.  Type  III  vehicles  were  not  included  in  the  data  analysis  because  their  exhaust 
emission  concentrations  were  not  measured  due  to  the  fact  that  the  emission  sources  were 
in  a  different  plane.  However,  their  number  and  flow  rate  were  monitored.  Figure  4.1.1 
summarizes  the  fleet  distribution  per  vehicle  class  and  Figure  4.1.2  shows  their  flow 
rates.  8,478  (50%)  were  type  I  vehicles,  4,829  (29%)  were  type  II  vehicles,  and  3,563 
(21%)  were  type  III  vehicles.  Overall  17.75  vehicles  passed  the  spectrometer  in  a 
minute.  The  flow  rates  were  8.92,  5.08  and  3.75  vehicles/min  for  type  I  vehicles,  type  II 
vehicles,  and  type  III  vehicles  respectively. 
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Figure  4.1.1.  Total  number  of  vehicles  monitored  during  PVEMP. 
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Figure  4.1.2.  Average  vehicle  flow  rate  recorded  during  PVEMP. 

The  measured  vehicle  speeds  were  integrated  in  order  to  calculate  the  speeds  and 
accelerations  at  the  spectrometer.  These  speeds  and  accelerations  were  also  correlated 
with  vehicle  type,  time  of  vehicle  at  the  spectrometer  and  time  of  exhaust  emission 
concentration  measurement.  Hourly  meteorological  data  were  then  added  to  the  other 
parameters  to  develop  the  PVEMP  database. 

The  integrated  final  PVEMP  database  consisted  of  4,413  vehicles.  2,848  (65%) 
were  type  I  vehicles  and  1,565  (35%)  were  type  II  vehicles  (Figure  4.1.3).  The  decrease 
in  the  number  of  vehicle  monitored  was  due  to  many  factors.  These  factors  included  the 
time  correlation  of  all  the  measured  parameters,  the  computational  time  required  for 
emission  concentration  calculation  by  the  FTIR  system,  the  missing  vehicle  clocking  for 
speed  measurement,  the  suspect  emission  concentration  values  that  were  discarded,  and 
the  subtraction  of  type  III  vehicles  from  the  analysis.  Appendix  A  shows  a  portion  and 
the  format  of  the  PVEMP  database. 


Figure  4.1.3.  Total  number  of  vehicles  included  in  the  data  analysis. 
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4.2.  Descriptive  Statistics 
4.2.1.  CO  and  HCs  Exhaust  Emission  Concentrations 

The  mean,  median,  minimum,  and  maximum  CO  and  HCs  exhaust  emission 
concentrations  for  all  the  vehicles  and  for  the  two  different  vehicle  types  are  presented  in 
figures  4.2.1  and  4.2.2. 
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Figure  4.2.1.  Mean  and  median  CO  and  HCs  exhaust  emission  concentrations  for  type  I 
and  type  II  vehicles,  and  the  combined  fleet. 
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Figure  4.2.2.  Minimum  and  maximum  CO  and  HCs  exhaust  emission  concentrations  for 
type  I  and  type  II  vehicles,  and  the  combined  fleet. 


The  mean  CO  exhaust  emission  concentrations  from  type  I  vehicles,  type  II 
vehicles,  and  the  combined  fleet  are  1.10  %,  1.16  %,  and  1.12  %  respectively.  The 
medians  are  0.63  %,  0.64  %  and  0.64  %  respectively.  This  suggests  that  the  distributions 
of  CO  exhaust  emission  concentration  are  severely  skewed  as  seen  in  Figure  4.2.3.  The 
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skewnesses  are  equal  to  4.05  %,  4.39  %,  and  4.21  %  for  type  I  vehicles,  type  II  vehicles, 
and  the  combined  fleet  respectively.  In  fact,  CO  concentration  displays  a  logarithmic 
distribution  as  depicted  in  Figures  4.2.4. 
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Figure  4.2.3.  Histogram  of  the  CO  exhaust  emission  concentration  distribution  for  the 
combined  fleet. 
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Figure  4.2.4.  Logarithmic  distribution  of  CO  exhaust  emission  concentration  for  the 
combined  vehicle  fleet  monitored. 


The  data  also  show  a  high  variability  in  CO  exhaust  emission  concentrations.  In  fact,  the 
minimum,  maximum,  and  standard  deviation  of  CO  exhaust  emission  concentration  for 
the  combined  fleet  are  0.01  %,  19.69  %,  and  1.61  %  respectively.  It  is  worth  mentioning 
that  Type  I  vehicles  registered  the  highest  CO  exhaust  emission  concentration  maximum 
and  the  lowest  mean  and  type  II  registered  the  lowest  CO  exhaust  emission  concentration 
minimum  and  the  highest  mean. 
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The  mean  HCs  exhaust  emission  concentrations  expressed  as  propane  equivalents 
from  type  I,  type  II,  and  the  combined  fleet  are  0.23  %,  0.24  %,  and  0.23  %  respectively. 
The  medians  are  0.19  %,  0.20  %  and  0.20  %  respectively.  This  suggests  that  the 
distributions  of  HCs  exhaust  emission  concentrations  are  somewhat  skewed  as  displayed 
in  Figure  4.2.5.  The  skewnesses  are  equal  to  1.28  %,  1.30  %,  and  1.30  %  for  type  I 
vehicles,  type  II  vehicles,  and  the  combined  fleet  respectively.  HCs  exhaust  emission 
concentration  also  displays  a  logarithmic  distribution  as  depicted  in  Figures  4.2.6.  The 
variability  of  the  HCs  exhaust  emission  concentration  distribution  is  lower.  The 
minimum,  maximum,  and  standard  deviation  of  HCs  exhaust  concentration  for  the 
combined  fleet  are  0.00  %,  1.40  %,  and  0.17  %  respectively.  Again,  type  I  vehicles 
registered  the  highest  HCs  exhaust  emission  concentration  maximum  while  type  II 
registered  the  highest  mean  HCs  exhaust  emission  concentration.  Appendix  B  presents 
the  histograms  and  logarithmic  distributions  for  each  vehicle  type. 

The  mean  concentrations  recorded  in  this  study  are  in  line  with  the  results  of  the 
few  remote-sensing  measurements  already  done.  For  example,  in  a  study  conducted  in 
Chicago,  Illinois  in  October  1991,  Stedman  and  co-workers  reported  mean  CO  and  HCs 
exhaust  emission  concentrations  of  1.10%  and  0.14%  respectively  (Stedman  et  al,  1991). 
In  the  study  conducted  in  El  Monte,  California  in  June  1991  by  General  Motor  Research 
and  Development  Center  and  the  California  Air  Resources  Board,  Stephens  reported 
mean  exhaust  emission  concentrations  of  0.86%  CO  and  0.04%  HCs  (Stephens,  1994). 
These  data  relate  to  the  combined  fleet  monitored  in  these  studies. 
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Figure  4.2.5.  Histogram  of  the  HCs  exhaust  emission  concentration  distribution  for  the 
combined  fleet. 
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Figure  4.2.6.  Logarithmic  distribution  of  HCs  exhaust  emission  concentration  for  the 
combined  vehicle  fleet  monitored. 


Figures  4.2.7  and  4.2.8  show  the  mean  CO  and  HCs  exhaust  emission 
concentrations  for  each  monitoring  session.  As  expected,  the  mean  concentrations  vary 
from  one  monitoring  session  to  the  other  due  to  changes  in  the  experimental  setups,  the 
driving  modes,  and  the  atmospheric  conditions.  The  highest  mean  CO  exhaust  emission 
concentrations  were  obtained  during  the  monitoring  session  of  May  20,  1997,  on  the 
Borman  Expressway.  The  lowest  were  obtained  during  that  of  August  21,  1997,  on  1-65. 
The  highest  mean  HCs  exhaust  emission  concentrations  were  obtained  during  the 
monitoring  session  of  May  20,  1997  and  the  lowest  during  that  of  August  22,  1997,  on  I- 
65.  The  physical  characteristics  of  these  sites  are  described  in  Section  3.  The 
atmospheric  conditions  on  these  days  are  described  below  in  the  present  Section. 

Figure  4.2.9  presents  the  percentiles  of  the  concentrations  of  CO  and  HCs  for  the 
combined  fleet.  The  plots  for  the  two  vehicle  types  are  presented  in  Appendix  B.  The 
interpretation  of  these  figures  is  critical  for  the  model  development  and  the  design  of 
concentration  threshold  limits  for  each  vehicle  type  because  they  can  report  on  the  data 
clusters  in  the  distributions.  For  example,  when  all  the  vehicles  are  considered,  it  can  be 
seen  that  90%  of  the  HCs  exhaust  emission  concentrations  are  less  than  0.46  %  and  10% 
are  less  than  0.04  %.  Furthermore,  90%  of  the  CO  exhaust  emission  concentrations  are 
less  than  2.22  %  and  10%  are  less  than  0.18  %.  This  type  of  finding  will  help  in  deciding 
which  data  to  include  in  the  development  of  the  exhaust  emission  concentration  model. 
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Figure  4.2.7.  Mean  CO  concentrations  for  type  I  vehicles,  type  II  vehicles,  and  the 
combined  fleet  recorded  during  each  monitoring  session. 
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Figure  4.2.8.  Mean  HCs  concentrations  for  type  I  vehicles,  type  II  vehicles,  and  the 
combined  fleet  recorded  during  each  monitoring  session. 
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Figure  4.2.9.  CO  and  HCs  exhaust  emission  concentration  percentiles  for  the  combined 
fleet. 


4.2.2.  Vehicle  Speed  and  Acceleration/Deceleration 

The  data  pertaining  to  the  driving  modes  is  summarized  below.  Vehicle  speeds 
varied  from  34.5  to  83.5  mph  with  a  standard  deviation  of  5.6  mph.  The  mean  vehicle 
speed  for  type  I  vehicles  was  56.74  mph  and  that  of  type  II  vehicles  was  56.21  mph.  It 
was  56.56  mph  considering  the  combined  fleet.  Figure  4.2.10  summarizes  the  data. 
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Figure  4.2.10.  Mean,  median,  maximum,  and  minimum  speeds  for  type  I  and  type  II 
vehicles,  and  the  combined  fleet. 
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The  vehicle  speed  distributions  exhibit  a  near  normal  characteristic  since  the  means  are 
almost  equal  to  the  medians.  The  day-to-day  variations  in  mean  vehicle  speeds  were 
between  55.2  mph  and  58.5  mph  depending  on  vehicle  type  (Figure  4.2.1 1). 
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Figure  4.2.1 1.  Mean  speeds  for  type  I  and  type  II  vehicles,  and  the  combined  fleet 
recorded  during  each  monitoring  session. 


During  the  monitoring  session  of  May  20,  1997,  when  the  highest  mean  CO  and  HCs 
exhaust  emission  concentrations  were  recorded,  the  mean  vehicle  speeds  for  type  I 
vehicles,  type  II  vehicles,  and  the  combined  fleet  were  56.9  mph,  57.3  mph,  and  57.0 
mph  respectively.  On  August  21,  1997,  when  the  lowest  mean  CO  exhaust  emission 
concentrations  were  recorded,  the  mean  vehicle  speeds  for  type  I  vehicles,  type  II 
vehicles,  and  the  combined  fleet  were  57.7  mph,  58.5  mph,  and  58.1  mph  respectively. 
On  August  22,  1997,  when  the  lowest  mean  HCs  exhaust  emission  concentrations  were 
recorded,  the  mean  vehicle  speeds  for  type  I  vehicles,  type  II  vehicles,  and  the  combined 
fleet  were  56.2  mph,  55.7  mph,  and  56.1  mph  respectively.  These  findings  are  also 
depicted  on  the  scatter  plots  of  the  mean  concentrations  versus  the  mean  vehicle  speeds 
for  the  combined  fleet  (Figure  4.2.12). 
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Figure  4.2.12.  Variation  of  mean  CO  and  HCs  exhaust  emission  concentration  with  mean 
vehicle  speed  for  the  combined  vehicle  fleet. 

The  acceleration  values  varied  from  -4.5  to  4.0  mph/sec  with  a  mean  of  0.17 
mph/sec,  a  median  of  0.0  mph/sec,  and  a  standard  deviation  of  0.71mph/sec  considering 
all  the  vehicles  (Figure  4.2.13). 
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Figure  4.2.13.  Mean,  median,  maximum,  and  minimum  acceleration  for  type  I  and  type  II 
vehicles  and  the  combined  fleet. 


The  day-to-day  variation  in  mean  acceleration  values  was  between  0.04  and  0.13 
mph/s  depending  on  car  type  except  on  August  21  and  22,  1997  when  the  mean 
acceleration  values  varied  from  0.36  to  0.59  mph/s  (Figure  4.2.14).  This  clearly  shows 
the  effect  of  the  characteristics  of  the  experimental  sites.  The  1-65  site  was  unobstructed 
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and  downhill  while  that  of  the  Borman  Expressway  was  a  flat  roadway  under  an 
overpass.  Except  for  August  22,  the  mean  accelerations  for  type  II  vehicles  were  higher 
than  those  of  type  I  vehicles.  Note  that  August  21  and  22,  1997  were  the  dates  when  the 
lowest  mean  CO  and  HCs  exhaust  emission  concentrations  were  respectively  recorded. 
The  variations  of  the  mean  concentrations  with  the  mean  vehicle  accelerations  are  also 
shown  in  Figure  4.2.15  for  the  combined  fleet. 
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Figure  4.2.14.  Mean  accelerations  for  type  I  and  type  II  vehicles,  and  the  combined  fleet 
recorded  during  each  monitoring  session. 
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Figure  4.2.15.  Variation  of  the  mean  CO  and  HCs  concentrations  with  the  mean  vehicle 
acceleration  for  the  combined  vehicle  fleet. 
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4.2.3.  Meteorology 

The  mean  meteorological  data  collected  during  the  field  experiments  is  presented 
in  Figure  4.2.16.  Ambient  temperature  varied  from  54.6  °F  to  72.0  °F  with  an  average  of 
61 .4  °F.  Wind  speed  varied  from  8.0  mph  to  12.5  mph  with  an  average  of  10.4  mph. 
Relative  humidity  varied  from  34.4%  to  79.9%  with  an  average  of  46.6%.  The  winds 
were  from  the  south-southeast  except  on  May  7,  1997  when  they  were  from  the 
northwest.  On  May  20,  1997,  when  the  highest  mean  CO  and  HCs  exhaust  emission 
concentrations  were  recorded,  the  mean  temperature,  wind  speed,  and  relative  humidity 
were  57.1  °F,  1 1.5  mph,  and  70.2%  respectively.  On  August  21,  1997,  when  the  lowest 
mean  CO  exhaust  emission  concentrations  were  recorded  they  were  62.6  °F,  10.0  mph, 
and  79.9%  respectively.  On  August  22,  1997,  when  the  lowest  mean  HCs  exhaust 
emission  concentrations  were  recorded,  they  were  72.0  °F,  12.4  mph,  and  50.0% 
respectively. 
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Figure  4.2.16.  Mean  meteorological  data  recorded  during  each  monitoring  session. 
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4.3.  Observed  Emission-Speed-Acceleration  Correlations 

The  interpretation  of  the  functional  relationships  between  real-time  CO  and  HCs 
exhaust  emission  concentrations,  vehicle  speeds,  driving  modes  (acceleration, 
deceleration,  and  cruise),  and  mean  values  of  atmospheric  parameters  is  very  important. 
It  is  the  key  to  the  development  of  a  sound  and  efficient  modal  exhaust  emission 
concentration  model  that  can  be  used  to  help  improve  traffic  flow  when  air  quality 
becomes  an  important  traffic  parameter.  The  first  approach  is  the  visual  study  of  the 
scatter  plots  of  the  key  parameters.  The  scatter  plots  of  the  mean  CO  and  HCs  emission 
concentrations  versus  the  mean  values  of  vehicle  speeds,  and  accelerations,  and 
atmospheric  observations  were  shown  in  Section  4.2.  The  scatter  plots  of  the 
instantaneous  values  of  CO  and  HCs  exhaust  emission  concentrations  versus  vehicle 
speeds  and  accelerations,  the  main  parameters  of  modal  emission  modeling,  are  presented 
in  Figures  4.3.1  to  4.3.4  for  the  combined  fleet.  Those  pertaining  to  the  different  vehicle 
types  are  presented  in  Appendix  B.  These  visual  correlations  do  not  show  any  classical 
mathematical  trends  except  the  fact  that  the  emission  concentration  data  are  clustered 
between  the  speeds  of  45.0  and  70.0  mph  and  the  accelerations  of -2.0  to  +  2.0  mph/sec. 
These  observations  may  be  due  to  many  factors.  Because  of  the  nature  of  the 
experimental  setups  and  the  characteristics  of  the  sites,  most  of  the  speed  and  driving 
mode  data  fall  in  the  above  speed  and  acceleration  ranges.  Moreover,  the  formation  of 
CO  and  unburned  HCs,  as  explained  in  Section  1,  makes  it  apparent  that  the  bulk  of  the 
exhaust  emissions  will  occur  in  these  ranges.  The  Pearson  correlation  factors  between 
the  instantaneous  values  of  CO  exhaust  emission  concentrations  and  vehicle  speeds  and 
accelerations  are  -0.046  and  0.005  respectively.  Those  between  the  instantaneous  values 
of  HCs  emission  concentrations  and  vehicle  speeds  and  accelerations  are  0.001  and  - 
0.033  respectively.  This  suggests  that  there  is  not  a  linear  relationship  between  emission 
concentrations,  speeds,  and  accelerations  and  neither  speed  nor  acceleration  alone  may 
explain  the  emission  concentration  variations.  Therefore  both  speed  and  acceleration 
should  be  accounted  for  simultaneously  when  studying  vehicle  exhaust  emission 
concentrations.  Figures  4.3.5  and  4.3.6  show  the  three-dimensional  plots  of  the  triple 
correlations  for  the  combined  fleet.  These  surfaces  were  derived  from  interpolated 
meshes  using  the  inverse  distance  method.  The  three-dimensional  plots  of  the  two 
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different  vehicle  types  are  presented  in  Appendix  B.  As  it  can  be  seen,  the  surfaces  are 
neither  uniform  nor  smooth  indicating  once  again  the  variability  and  cluster  of  the 
exhaust  emission  concentrations.  The  determination  of  these  three-dimensional 
functional  relationships  will  represent  the  theoretical  basis  of  the  model. 
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Figure  4.3.1.  Variation  of  instantaneous  CO  exhaust  emission  concentrations  with  vehicle 
speed  for  the  combined  vehicle  fleet. 
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Figure  4.3.2.  Variation  of  instantaneous  CO  exhaust  emission  concentrations  with  vehicle 
acceleration  for  the  combined  vehicle  fleet. 
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Figure  4.3.3.  Variation  of  instantaneous  HCs  exhaust  emission  concentrations  with 
vehicle  speed  for  the  combined  vehicle  fleet. 
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Figure  4.3.4.  Variation  of  instantaneous  HCs  exhaust  emission  concentrations  with 
vehicle  acceleration  for  the  combined  vehicle  fleet. 
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Figure  4.3.5.  Three-dimensional  plot  for  the  triple  correlation  of  CO  exhaust  emission 
concentration,  speed,  and  acceleration  for  the  combined  fleet. 
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Figure  4.3.6.  Three-dimensional  plot  for  the  triple  correlation  of  HCs  exhaust  emission 
concentration,  speed,  and  acceleration  for  the  combined  fleet. 


4.4.  High  Emitters 

The  CO  and  HCs  high  emitting  vehicle  fractions  along  with  their  contribution  to 
the  total  exhaust  emission  concentrations  are  depicted  in  Figures  4.4.1  and  4.4.2 
respectively  for  the  combined  fleet.  The  cut-points  for  high  emitters  were  selected  to  be 
4.0  %  CO  and  0.2  %  HCs  (Jack  et  al,  1995,  Stephens,  1994).  The  fractions  of  the  total 
vehicles  emitting  at  or  above  a  given  emission  concentration  level  are  also  presented  in 
Figures  4.1.3  and  4.1.4.  Overall,  5.7%  of  the  vehicle  fleet  emitted  more  than  4%  CO  and 
generated  33.7%  of  the  total  CO  exhaust  emission  concentration.  49.2%  of  the  vehicle 
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fleet  emitted  more  than  0.2%  HCs  and  generated  77.5%  of  the  total  HCs  exhaust 
emission  concentration.  However,  Only  3.0%  of  the  total  vehicle  fleet  were  high  emitters 
of  both  CO  and  HCs  together. 
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Figure  4.4.1.  Fractions  of  CO  high  emitters  and  total  CO  exhaust  emission  concentrations 
for  the  combined  vehicle  fleet. 


0.80 


E]Veh.  Fraction 
■  HC  Fraction 


Type  I  Type  II  All 

Vehicle  Type 


Figure  4.4.2.  Fractions  of  HCs  high  emitters  and  total  HCs  exhaust  emission 
concentrations  for  the  combined  vehicle  fleet. 
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5.  REAL-TIME  MODAL  MODEL  DEVELOPMENT 

5.1.  Modeling  Overview 

Models  are  formal  set  of  calculations  designed  to  closely  represent  reality.  Their 
main  roles  are  to  predict  or  describe  unmeasured  conditions,  predict  and  examine  changes 
and  causalities,  and  avoid  the  high  costs  of  data  collection  and  monitoring.  The  major 
components  of  a  model  are  the  theories,  the  input  parameters  and  data,  and  the  underlying 
assumptions.  Since  models  are  set  of  calculations,  they  bear  inherent  errors  rising  from 
the  assumptions,  the  limitations  of  the  theories,  and  the  input  data. 

When  modeling  traffic  emissions,  three  popular  computational  methods  are  used. 
They  are  the  method  based  on  real  driving  behavior,  the  method  based  on  road  stretch 
analysis,  and  the  method  based  on  mileage  related  emission  balances.  The  scale  of 
applications,  the  levels  of  details,  and  the  input  data  requirements  of  these  computational 
method  are  different.  The  emission  calculations  based  on  actual  driving  behavior  (modal 
modeling)  are  applied  on  a  local  (micro)  scale.  The  emission  factors  are  emission  maps 
for  vehicles  based  on  modal  data  from  transient  test  or  real  time  monitoring.  The  traffic 
data  are  derived  from  counts  or  traffic  models.  The  road  specific  data  are  represented  by 
actual  driving  behaviors.  The  emission  calculations  based  on  road  stretch  analysis  are 
applied  on  a  regional  scale.  The  emission  factors  are  generally  from  FTP  dynamometer 
tests.  The  traffic  data  are  derived  from  counts  or  traffic  models.  Mean  speed,  road 
gradient,  and  altitude  for  each  road  represent  the  road  specific  data.  The  emission 
calculations  based  on  mileage-related  emission  balance  are  applied  on  a  regional, 
national,  and  global  scale.  The  emission  factors  are  generally  from  FTP  dynamometer 
tests.  The  traffic  data  are  derived  from  statistical  data  for  VMTs  according  to  subclasses 
of  vehicles  and  roads  and  there  is  no  road  specific  data. 

Modal  emission  modeling  is  relatively  new  and  encompasses  three  different 
methods.  A  convenient  method  to  characterize  vehicle-operating  modes  is  the 
development  of  speed-acceleration  matrix  that  measures  emissions  associated  with  each 
mode.  Another  modal  emissions  method  develops  an  emission  map  based  on  engine 
power  and  speed.  These  two  methods  have  used  second-by-second  emission  tests 
performed  at  numerous  engine-operating  points,  taking  an  average  of  steady-state 
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measurements.  St-denis  and  Winer  (1993)  have  created  both  a  speed-acceleration  and  a 
speed-load  modal  emissions  model  using  data  from  a  single  Ford  vehicle.  Sierra 
Research  has  developed  a  modal  emissions  model  that  computes  the  second-by-second 
engine  speed  and  load  required  for  a  specified  driving  cycle,  then,  using  an  emissions 
map,  second-by-second  emissions  were  approximated  (Barth  et  al,  1998).  The  other 
modal  emission  method  is  the  analytical  power-demand  modeling.  It  is  based  on 
parameterized  analytical  representation  of  emission  production.  In  this  extensive 
treatment,  the  entire  emission  process  is  broken  down  into  components  that  correspond  to 
physical  phenomena  associated  with  vehicle  operation  and  emissions  production.  Each 
component  is  then  modeled  as  an  analytical  representation  consisting  of  various 
parameters  that  are  characteristic  of  the  process.  These  parameters  vary  according  to 
vehicle  type,  engine,  and  emission  technology.  Researchers  at  the  University  of 
California,  Riverside,  are  currently  developing  such  a  model  (Barth  et  al,  1998). 
However,  so  far,  on-road  real-time  data  have  not  yet  been  used  in  modal  emissions 
modeling.  Moreover,  modal  emission  modeling  has  not  yet  been  used  for  the  estimation 
of  real-time  emission  production  based  on  input  parameters  from  real-time  advanced 
traffic  monitoring  systems. 

5.2.  Model  Development 
5.2.1.  Description  of  the  current  model 

The  model  developed  in  this  study;  the  Purdue  Vehicle  Exhaust  Emission  Model, 
the  Borman  Expressway  Application  (PVEM-BEA)  is  a  real-time,  micro-scale,  modal 
exhaust  emission  concentration  model.  It  is  designed  to  predict  average  vehicular 
exhaust  emission  concentrations  of  CO  and  HCs  in  real-time  traffic  flow.  It  is  a  data 
driven  model  based  on  real  driving  behavior  and  derived  from  the  integrated  CO  and  HCs 
exhaust  emission  concentrations  in  terms  of  the  modal  value-pairs  speeds  and 
accelerations.  The  regression  analyses  of  the  speed-acceleration  matrices  result  from  the 
modal  analysis  of  the  PVEMP  database  described  in  Section  4.  PVEM-BEA  is  intended 
to  be  an  integral  part  of  the  Borman  Expressway  Advanced  Traffic  Management  System 
representing  an  additional  tool  for  traffic  management,  especially  in  construction  zones, 
when  air  pollution  becomes  an  important  decision-making  parameter. 
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Real-time  values  of  vehicle  speed  and  acceleration,  traffic  composition  and 
volume,  and  road  geometry  are  the  main  input  parameters.  They  are  the  main  parameters 
affecting  real-time  vehicle  exhaust  emissions.  Even  though  vehicle  age,  conditions,  and 
fuel  economy  affect  vehicle  exhaust  emissions  also,  they  are  excluded  in  PVEM-BEA 
because  they  can  not  be  assessed  instantaneously.  Synoptic  meteorology  is  also  excluded 
because  of  the  small  spatial  and  temporal  scale  of  the  model. 

For  a  given  location,  the  real-time  input  of  speed,  acceleration,  and  type  of  each 
passing  vehicle  acquired  instantaneously  (very  short  time  interval)  are  used  to  calculate 
the  mode-average  CO  and  HCs  exhaust  emission  concentrations  according  to  the  vehicle 
type.  The  program  then  outputs  data  about  the  traffic  volume,  the  fleet  composition,  the 
vehicle  flow  rates,  the  vehicle  speeds  and  accelerations,  the  average  CO  and  HCs  exhaust 
emission  concentrations,  and  the  emission  concentration  threshold  level  exceedences 
(decision  making)  for  the  short  time  interval. 

5.2.2.  Analytical  Background  of  PVEM-BEA 

It  was  shown  in  Section  4  that  there  were  no  mathematical  relationships  between 
the  instantaneous  CO  and  HCs  exhaust  emission  concentration  and  speeds  and 
accelerations.  Consequently,  the  database  has  to  be  reduced  further  in  order  to  assess  the 
correlations.  In  doing  so,  speeds  and  accelerations  were  binned  and  average  exhaust  CO 
and  HCs  emission  concentrations  were  calculated  at  acceleration  midpoints  independent 
of  speed,  at  speed  midpoints  independent  of  acceleration,  and  at  the  intersection  of  the 
speed-acceleration  pairs.  The  acceleration  bins  were  constructed  as  follows:  the 
midpoints  were  selected  to  be  -4.0,  -3.5,  -3.0,  ...,  -0.5,  0.0,  0.5,  ...,3.0,  3.5, 4.0  mph/sec. 
The  ranges  were  selected  to  be  [(midpoint  -  0.25);  (midpoint  +  0.25)[  for  acceleration  and 
](midpoint-0.25);  (midpoint  +  0.25)]  for  deceleration.  Speed  values  were  varied  at 
increment  of  1.0  mph  from  30.0  mph  to  80.0  mph.  This  methodology  was  based  on  the 
resolution  of  the  measured  speed  values.  Figures  5.2.1  to  5.2.4  show  the  plots  of  the  bin- 
average  CO  and  HCs  exhaust  emission  concentrations  versus  speed  and  acceleration 
independent  of  each  other  for  the  combined  fleet.  The  same  plots  were  developed  for  the 
two  vehicle  types  and  are  presented  in  Appendix  C.  Both  bin-average  CO  and  HCs 
exhaust  emission  concentrations  exhibit  high  variabilities  at  both  ends  of  the  speed  and 
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acceleration  ranges;  that  are  high  acceleration  and  deceleration,  and  high  and  low  speeds. 
Within  the  common  driving  mode  ranges  of  the  Borman  Expressway,  i.e.,  speeds  within 
45.0  to  70.0  mph  and  accelerations  within  -2.0  and  2.0  mph/sec,  the  bin-average  CO  and 
HCs  exhaust  emission  concentrations  are  nearly  linear. 
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Figure  5.2.2.  Bin-average  HCs  exhaust  emission  concentration  versus  speed  midpoints. 
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Tigure  5.2.3.  Bin-average  CO  exhaust  emission  concentration  versus  acceleration 
midpoints. 
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Figure  5.2.4.  Bin-average  HCs  exhaust  emission  concentration  versus  acceleration 
midpoints. 
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Second  order  polynomials  were  fit  to  the  data  as  shown  in  the  Figures.  Because  of  the 
high  variabilities  at  both  ends  of  the  speed  and  acceleration  spectra,  the  coefficients  of 
correlation  were  very  low  ranging  from  20  to  40%.  This  suggests  that  neither  speed  nor 
acceleration  alone  may  explain  all  the  variations  of  the  bin-average  CO  and  HCs  exhaust 
emission  concentrations.  However,  in  real-time  monitoring,  only  these  two  variables 
may  be  assessed.  Consequently,  the  two  variables  should  be  used  together  when 
analyzing  the  bin-average  CO  and  HCs  exhaust  emission  concentrations. 

Using  the  acceleration  and  speed  bins  developed,  the  average  CO  and  HCs 
exhaust  emission  concentrations  were  calculated  at  each  speed-acceleration  value  pairs 
(or  modes).  This  defined  a  speed-acceleration  matrix  for  mode-average  CO  and  HCs 
exhaust  emission  concentrations  (Figure  5.2.5).  Figures  5.2.6  and  5.2.7  show  the  three- 
dimensional  plots  of  the  mode-average  CO  and  HCs  exhaust  emission  concentrations 
against  speed  and  acceleration  midpoints.  The  surfaces  were  also  derived  from 
interpolated  meshes  using  the  inverse  distance  method.  The  three-dimensional  plots  of 
the  two  different  vehicle  types  are  presented  in  Appendix  C.  As  it  can  be  seen,  the 
surfaces  are  neither  uniform  nor  smooth  indicating  that  the  mode-average  CO  and  HCs 
exhaust  emission  concentrations  are  also  variable  and  clustered. 
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Figure  5.2.5.  Speed-acceleration  matrix  for  Eij  (Eij  =  mode-average  exhaust  emission 
concentration  at  speed  i  and  acceleration  j). 
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Figure  5.2.6.  Three-dimensional  plot  of  the  mode-average  exhaust  CO  emission 
concentration  versus  speed  and  acceleration  midpoints  for  the  combined  fleet. 
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Figure  5.2.7.  Three-dimensional  plot  of  the  mode-average  exhaust  HCs  emission 
concentration  versus  speed  and  acceleration  midpoints  for  the  combined  fleet. 


Using  a  non-linear  mutivariate  regressive  approach,  the  mode-average  CO  and 
HCs  exhaust  emission  concentrations  variations  with  speed  and  acceleration  midpoints 
could  be  smoothed  out  for  estimation  purposes.  Based  on  the  variations  of  Figures  5.2.1 
to  5.2.4,  a  second  order  multiple  regression  on  both  speed  and  acceleration  was  used. 
The  interactions  between  speed  and  acceleration  were  also  included.  The  Marquardt- 
Levenberg  algorithm  was  used  to  find  the  coefficients  (parameters)  of  the  independent 
variables  (speed  and  acceleration)  that  give  the  best  fit  between  the  equation  and  the  data. 
This  algorithm  seeks  the  values  of  the  parameters  that  minimize  the  sum  of  the  squared 
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differences  between  the  values  of  the  observed  and  predicted  values  of  the  dependent 
variable.  This  process  is  iterative.  It  begins  with  a  guess  at  the  parameters,  checks  to  see 
how  well  the  equation  fits,  then  continues  to  make  better  guesses  until  the  differences 
between  the  residual  sum  of  squares  no  longer  decreases  significantly  or  converges.  The 
following  predictive  equation  for  the  mode-average  CO  and  HCs  exhaust  emission 
concentrations  was  determined  based  on  this  method. 

E  =  CO  +  Cl*((S-55.0)/30.0)  +  C2*((S-55.0)/30.0)2 
+  C3*(A/4.0)  +  C4*(A/4.0)2  +  C5*((S-55.0)/30.)*(A/4.0)  (1) 

where  CO,  CI,  C2,  C3,  C4,  and  C5  are  the  regression  coefficients.  E  represents  the 
mode-average  CO  and  HCs  exhaust  emission  concentrations  in  %,  S  the  speed  in  mph, 
and  A  the  acceleration  in  mph/sec.  Tables  5.2.1  and  5.2.2  show  the  values  of  the 
regression  coefficients  and  that  of  the  multiple  coefficient  of  determination  for  the  mode- 
verage  CO  and  HCs  exhaust  emission  concentrations  for  the  different  vehicle  classes. 
The  multiple  coefficient  of  determination  is  defined  as 

R2  =  sum[(0-P)2]/sum[(0-avg(0))2]  (2) 

O  and  P  are  the  observed  and  predicted  mode-verage  CO  or  HCs  exhaust  emission 
concentrations.  The  very  high  values  of  the  multiple  coefficient  of  determination  implies 
that  the  model  fits  very  well  the  data  and  that  at  least  96%  of  the  variation  in  the  mode- 
average  CO  and  HCs  exhaust  emission  concentrations  are  attributable  to  speed  and 
acceleration. 


Vehicle  types 

Parameters 

Type  I 

Type  II 

All 

CO 

1.2490 

1.3240 

1.1780 

CI 

-0.2855 

0.0908 

-0.2403 

C2 

-0.6823 

-1.0890 

-0.6993 

C3 

0.0013 

0.1417 

0.0225 

C4 

0.2185 

-0.2069 

0.4413 

C5 

1.0440 

-0.0418 

0.4337 

R2 

0.9792 

0.9934 

0.9682 

Table  5.2.1.  Regression  parameters  for  the  mode-average  CO  exhaust  emission 
concentration  against  speed  and  acceleration. 
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Vehicle  types 

Parameters 

Type  I 

Type  II 

All 

CO 

0.2324 

0.2471 

0.2367 

CI 

-0.0231 

-0.0385 

-0.0261 

C2 

0.0080 

-0.0494 

-0.0068 

C3 

-0.0274 

-0.0396 

-0.0278 

C4 

0.0373 

0.0233 

0.0292 

C5 

0.0539 

-0.0375 

-0.0330 

R* 

0.9905 

0.9799 

0.9841 

Table  5.2.1.  Regression  parameters  for  the  mode- average  HCs  exhaust  emission 
concentration  against  speed  and  acceleration. 


Figures  5.2.8  and  5.2.9  show  the  three-dimensional  plots  of  the  application  of  the 
mathematical  relationships  for  the  combined  fleet.  Those  for  the  two  different  vehicle 
types  are  included  in  Appendix  C.  These  surfaces  were  also  derived  from  interpolated 
meshes  using  the  inverse  distance  method.  As  it  can  be  seen,  the  mode-average  CO  and 
HCs  exhaust  emission  concentrations  are  highest  at  hard  accelerations  and  decelerations 
and  flatten  off  at  mild  accelerations.  This  finding  is  in  agreement  with  the  theories  of  the 
formation  of  CO  and  HCs  exhaust  emissions  described  in  the  Introduction  section  (Turns, 
1996).  It  can  also  be  seen  that  the  CO  surfaces  are  smoother  than  the  HCs  surfaces. 
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Figure  5.2.8.  Three-dimensional  plot  of  the  mode-average  CO  exhaust  emission 
concentration  obtained  from  the  mathematical  relationship  for  the  combined  fleet. 
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Figure  5.2.9.  Three-dimensional  plot  of  the  mode-average  HCs  exhaust  emission 
concentration  obtained  from  the  mathematical  relationship  for  the  combined  fleet. 


5.2.3.  Limitations  and  Extrapolations  of  PVEM-BEA 

The  mathematical  relationships  built  from  the  multivariate  regression  approach 
may  not  be  applied  outside  the  ranges  of  the  supporting  experimental  data  that  are  speeds 
between  30.0  mph  and  80.0  mph  and  accelerations  between  -4.0  mph/sec  and  4.0 
mph/sec.  Since  the  model  is  primarily  intended  to  be  used  at  construction  zones, 
extrapolation  of  the  estimation  method  is  necessary  in  order  to  include  lower  speed 
values.  The  extrapolation  method  used  in  this  study  depends  solely  on  acceleration  since 
it  represents  the  critical  and  limiting  variable  as  far  as  real-time  estimation  of  exhaust 
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emission  concentrations  is  concerned.  Therefore,  for  speed  values  outside  the 
experimental  data  range,  the  model  will  switch  to  an  acceleration  alone  mode.  Thus  the 
best  fit  of  the  mode-average  CO  and  HCs  exhaust  emission  concentrations  against 
acceleration  alone  (Figures  5.2.3  and  5.2.4)  given  in  Equation  2  will  be  used  for  the 
estimations 

E  =  CO  +  CI  *(A/4.0)  +  C2*(A/4.0)2  (3) 

Tables  5.2.3  and  5.2.4  summarize  the  parameters  of  Equation  3. 


Vehicle  types 

Parameters 

Type  I 

Type  II 

All 

CO 

1.1890 

1.173 

1.0830 

CI 

-0.3745 

-0.2512 

-0.2112 

C2 

0.5304 

0.4121 

0.2590 

Table  5.2.2.  Regression  parameters  for  the  mode-average  CO  exhaust  emission 
concentration  against  acceleration  alone. 

Vehicle  types 

Parameters 

Type  I 

Type  II 

All 

CO 

0.2357 

0.2293 

0.2315 

CI 

-0.0243 

0.0020 

-0.0124 

C2 

0.0703 

0.1026 

0.0766 

Table  5.2.3.  Regression  parameters  for  the  mode-average  HCs  exhaust  emission 
concentration  against  acceleration  alone. 


5.2.4.  Design  of  Concentration  thresholds  for  decision-making 

For  individual  vehicles,  the  design  of  the  CO  and  HCs  concentration  thresholds 
was  based  on  the  statistical  characteristics  of  the  regression  analysis  of  the  triple 
correlation  between  exhaust  emission  concentrations,  speeds,  and  accelerations.  The 
instantaneous  values  of  the  mode-average  CO  and  HCs  exhaust  emission  concentration 
thresholds  were  set  to  the  upper  limits  of  the  99%  confidence  intervals  of  the  regression 
means.  Thus,  for  individual  type  I  vehicles,  the  exhaust  emission  concentration 
thresholds  were  1.24  %  CO  and  0.24  %  HCs.  For  individual  type  II  vehicles,  they  were 
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1.26  %  CO  and  0.25  %  HCs.  For  the  combined  fleet  they  were  1.25  %  CO  and  0.24  % 
HCs. 

For  classes  of  vehicles  modeled  within  the  maximum  averaging  time,  the  design 
method  was  different.  Here,  the  measured  99th  percentile  of  CO  and  HCs  exhaust 
emission  concentrations  were  multiplied  by  the  measured  vehicle  flow  rates.  This 
product  was  called  the  concentration  threshold  limit  (CTL)  and  expressed  in  %/sec.  The 
product  of  the  estimated  class-average  CO  and  HCs  exhaust  emission  concentrations  by 
the  flow  rate  and  the  total  number  of  traffic  lanes  are  compared  to  the  CTLs  in  the  model 
to  arrive  at  a  decision  about  concentration  levels  within  the  maximum  averaging  time. 
For  type  I  vehicles,  the  CTLs  were  1.16  %/sec  CO  and  0.1 1  %/sec  HCs.  For  type  II 
vehicles,  they  were  0.68  %/sec  CO  and  0.06  %/sec  HCs.  For  the  combined  fleet  they 
were  1.86  %/sec  CO  and  0.17  %/sec  HCs. 

5.2.5.  Computer  Algorithm  of  PVEM-BEA 

The  current  version  of  PVEM-BEA  is  written  in  FORTRAN  for  the  PC 
computing  environment.  It  requires  a  FORTRAN  compiler  such  as  the  Microsoft 
FORTRAN  Visual  Workbench™  for  its  execution.  It  can  be  run  as  a  stand-alone 
program  or  as  integrated  in  the  Borman  Expressway  Traffic  Management  Center  (TMC). 
It  also  requires  input  from  Advanced  Traffic  Monitoring  systems  such  as  the 
Autoscope™  or  two  consecutive  pre-installed  conventional  inductive  loops.  Figure  5.2.10 
shows  a  top  down  flow  chart  of  PVEM-BEA  execution  algorithm.  Figure  5.2.1 1  shows  a 
sample  output  of  PVEM-BEA  obtained  from  a  hypothetical  simulation  input.  The  main 
program  makes  a  sequential  call  often  different  subroutines  whose  operations  are 
described  below. 

INCLUDE:  this  not  a  subroutine  but  a  call  to  an  included  file  that  initializes  all  the 
constant  parameters  including  the  regression  parameters  and  CO  and  HCs  emission 
concentration  thresholds. 
SUBROUTINES: 

OPENER:  Opens  all  the  input  and  output  files. 

READER1:  Inputs  the  simulation  location  name,  the  number  of  traffic  lanes,  and  the 
distance  between  loops  or  consecutive  speed  measurement  points. 
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WRITER1 :  Outputs  the  model  description,  and  the  variables  read  by  READER1 . 

READER2:  Inputs  real-time  traffic  parameters  including  speeds,  and  vehicle  type. 

Calculates  traffic  volumes  and  flow  rates  by  vehicle  class  within  the  maximum  time 

interval. 

CALCACCEL:  Calculates  average  speed  and  acceleration  of  each  passing  vehicle  within 

the  time  interval  of  speed-readings,  and  average  speed  and  acceleration  by  vehicle  type 

within  the  maximum  time  interval. 

CALCEMISSION:  Calculates  the  mode-average  CO  and  HCs  exhaust  emission 

concentration  for  each  passing  vehicle  as  well  as  for  classes  of  vehicles  within  the 

maximum  time  interval.  Compares  instantaneous  emission  concentrations  for  each 

vehicle  to  threshold  values.  Calculates  the  product  of  flow  rates  and  average  emission 

concentration  per  vehicle  type  for  decision-making  purposes. 

DECMAKER:  Makes  decision  about  average  CO  and  HCs  exhaust  emission 

concentrations  by  comparing  the  product  of  flow  rates  and  average  concentration  per 

vehicle  type  to  designed  threshold  levels  of  the  same  product. 

WRITER2:  Outputs  all  the  relevant  information  concerning  calculated  traffic  parameters, 

exhaust  emission  concentration,  decision  made  within  the  maximum  time  interval  for 

each  vehicles  as  well  as  for  each  class  of  vehicles. 

RESETER:  Loops  through  the  program  for  the  nest  simulation  time  or  calls  an  end  to  the 

program  if  the  simulation  is  over. 

CLOSER:  Closes  all  the  input  and  output  files  when  the  simulation  ends. 
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Figure  5.2.10.  Top-down  flow  chart  of  the  PVEM-BEA  subroutines. 
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********************************************************************************* 

*  Purdue  Real-Time  Vehicle  Exhaust  Emission  Modal  Model  V  1 .0  (PVEM1 .0) 

*  Developed  by  Dr.  Ouattara  Fatogoma.  Civil  Engineering,  Purdue  University.  All  Rights  Reserved. 
♦♦♦it*************************************************************************** 

NOTES 


CO  and  HCs  instantaneous  exhaust 
emission  concentration  thresholds(%) 
VT    CO       HCs 

1  1.3300  .3300 

2  1.3500  .3300 

3  .0000  .0000 

4  1.3400  .3300 

Designed  CO  and  HCs  threshold  limit  values 
by  vehicle  type(%/sec) 
VT    CO      HCs 

1  1.2500   .2000 

2  .7700  .1500 

3  .0000  .0000 

4  1.9500   .2600 

Notes:  bt  =  below  threshold,  at  =  above  threshold 

l=type  1  vehicles  (LDV),  2=type  2  vehicles  (MDV),  3=type  3  vehicles  (HDV),  4=Combined  1&2 

PCFNL=product  concentration,  flow  rate,  number  of  lanes 

Modeling  location:  Borman  Expressway  Mileage  100 

Number  of  lanes:  3 

Distance  between  2  speed  detection  points  (ft):  100.00 


Run#:  000001 

SIMULATION  DATA 

Date  Time  Vt   Spd    Ace 

(mmddyy)  (Thhrnmss)        (mph)  (mph/s) 

020998    '  120000  3    50.00      .00 

020998  120000  2    57.50   -1.50 

020998  120000  1    57.50    -4.00 

SUMMARY  DATA 

Number  of  vehicles: 
Vehicle  Flow  rate(#veh/sec/lane): 
Average  speed(mph): 
Average  acceleration(mph/sec): 
Average  CO  emission  concentration(%): 
Average  HC  emission  concentration(%): 
Number  of  vehicles  exceeding  CO  threshold: 
Number  of  vehicles  exceeding  HC  threshold: 
CO  PCFNL(%/sec): 
HCs  PCFNL(%/sec): 
CO  PCFNL  flag: 
HCs  PCFNL  flag: 


CO      COfg 

HC      HCfg 

(%) 

(%) 

.0000    bt 

.0000    bt 

1.2431     bt 

.2628    bt 

1.3507    at 

.2907   bt 

Type  1 

Type  2 

Type  3 

CF  1&2 

1 

1 

1 

2 

.08 

.08 

.08 

.17 

57.50 

57.50 

50.00 

57.50 

^1.00 

-1.50 

.00 

-2.75 

1.3507 

1.2431 

.0000 

1.2969 

.2907 

.2628 

.0000 

.2768 

:          1 

1 

0 

0 

:          1 

1 

0 

0 

.3377 

.3108 

.0000 

.6484 

.0727 

.0657 

.0000 

.1384 

bt 

bt 

bt 

bt 

bt 

bt 

bt 

bt 

Figure  5.2.11.  Sample  output  of  PVEM-BEA  for  single  simulation 
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6.  IMPLEMENTATION  SUGGESTIONS 

The  PVEM-BEA  Model  is  used  to  estimate  the  exhaust  emission  concentrations 
of  CO  and  HCs  in  real-time  traffic  flow.  The  only  inputs  required  of  the  model  are  two 
consecutive  speed-readings  for  a  single  vehicle  and  the  vehicle  type  identification.  The 
output  from  the  model  is  the  input  traffic  parameters,  the  emissions  concentrations  with 
decision  making  threshold  trigger,  and  the  time-average  data  for  each  class  of  vehicle. 

With  these  requirements  in  mind,  the  model  can  be  implemented  in  a  variety  of 
ways.  It  can  be  used  prior  to  a  project  to  assess  the  possible  CO  and  HCs  impacts  of  a 
construction  project  lane  closure,  during  a  project  to  confirm  the  impact  and  to  monitor 
the  situation,  or  after  the  fact  as  an  emissions  inventory  or  airshed  impact  tool.  Data  from 
current  Advanced  Traffic  Management  System  (ATMS)  can  be  transmitted  to  a  central 
station  where  they  can  be  processed,  or  the  data  can  be  collected  in-situ  and  processed 
real-time  via  a  laptop  computer. 

The  ATMS  currently  under  development  on  the  Borman  Expressway  can 
seamlessly  integrate  the  PVEM-BEA  model  into  the  network.  Any  data  that  is  collected 
and  stored  as  a  database  of  specified  format  can  be  read  into  the  model  to  produce  an 
emissions  concentration  estimate.  It  may  also  be  possible  to  read  the  data  into  the  model 
real-time  once  the  system  is  up  and  running. 

In  the  field,  the  Hoosier  Helper  vehicles  equipped  with  a  video  camera  on  a 
vertically  telescoping  mount  can  be  implemented  in  conjunction  with  an  Autoscope™ 
system  made  by  Image  Sensing  Systems  Incorporated.  The  method  of  detection  used  by 
this  system  requires  no  permanent  detector  installation  and  no  disruption  of  traffic  or  the 
roadway.  This  allows  the  system  to  be  highly  mobile  and  be  transported  to  the  site  of 
concern  for  quick  and  easy  determination  of  traffic  parameters.  Since  construction  zones 
are  of  main  concern,  this  mobility  is  of  critical  importance. 

With  the  Autoscope™  system,  a  video  image  is  analyzed  up  to  30  times  per  second 
and  is  capable  of  determining  vehicle  speed,  count,  and  type.  The  Autoscope™  unit 
shown  in  Figure  6.1,  detects  vehicle  presence  by  monitoring  the  video  image  at  a 
specified  location  on  the  roadway.  These  locations  are  essentially  sensors  on  the  roadway 
defined  by  the  interface  computer.  The  video  image  of  the  roadway  is  displayed  on  the 
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interface  computer  screen  and  boxes  are  drawn  with  the  mouse  to  identify  regions  where 
the  image  is  to  be  analyzed  and  then  identified  for  the  type  of  analysis  to  be 
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Figure  6.1.  The  Autoscope™  Video  Detection  Processor  Unit. 

conducted  (see  Figure  6.2).  Changes  in  the  image  attributes  at  these  specified  locations 
are  used  to  identify  and  measure  the  vehicle  type,  count,  and  speed.  The  detector  regions 
will  indicate  a  detected  vehicle  visually  on  the  computer  screen  by  changing  color  during 
setup  to  assure  proper  operation  (see  Figure  6.2). 


Vehicle  Type&Count  Detector 


First  Speed  Detector 


Second  Speed  Detector 
/ 


Figure  6.2.  Schematic  of  Video  Detection. 


Data  can  be  collected  for  multiple  lanes  on  either  an  instantaneous  or  a  time- 
averaged  basis.  The  instantaneous  data  can  be  read  directly  into  the  model  for 
instantaneous  vehicle  emissions  determination.  The  interval  data  is  a  time  average 
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dataset  and  can  be  generated  for  time  intervals  of  10,  20,  and  30sec,  or  1,  5,  10,  15,  30, 
and  60  minutes.  In  all  cases,  the  acceleration  parameter  will  be  calculated  within  the 
PVEM-BEA  model  from  multiple  speed-readings  of  the  same  vehicle.  This  in-situ 
method  of  data  collection  can  give  a  real-time  evaluation  of  the  effects  of  both  local 
ATMS  efforts  and  any  traffic  obstruction  such  as  an  incident  or  a  construction  zone. 
The  PVEM-BEA  model  can  be  a  valuable  planning  tool  as  well.  Prior  to  the 
beginning  of  a  construction  project,  vehicle  data  can  be  gathered  and  modified  to 
represent  the  event  of  a  lane  closure.  Once  modified,  the  dataset  can  be  input  into  the 
model  to  estimate  the  emissions  impact  of  the  lane  closure  or  speed  reduction.  In  this 
way,  multiple  scenarios  can  be  analyzed  and  decisions  made  to  minimize  the  actual 
impact  of  the  traffic  pattern  change. 
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7.  CONCLUSION 

The  Midac  FTIR  air  monitoring  system  was  used  to  measure  real-time  on-road 
vehicular  CO  and  HCs  exhaust  emission  concentrations  in  spring  and  summer  1997. 
Concurrently,  two  vehicle  speed-readings  were  acquired  from  the  passing  vehicles  in  the 
vicinity  of  the  spectrometer  using  a  hand-held  laser  gun.  The  speed  data  were  integrated 
to  calculate  average  vehicle  speed  and  acceleration  at  the  spectrometer.  Furthermore,  an 
8-mm  camcorder  was  set  up  to  videotape  the  passing  vehicles.  These  video  imageries 
were  reduced  in  order  to  classify  the  passing  vehicles.  Meteorology  data  were  also 
acquired  from  the  nearest  surface  meteorological  stations. 

Overall,  16,870  vehicles  were  monitored.  8,478  were  type  I  vehicles 
(automobiles),  4,829  were  type  II  (medium-duty  vehicles),  and  3,563  were  type  III 
(heavy-duty  vehicles).  The  flow  rates  of  the  vehicles  in  count/min,  were  8.92  for  type  I 
vehicles,  5.08  for  type  II,  3.75  for  type  HI,  and  17.75  for  the  combined  fleet.  4,413 
vehicles  were  retained  for  data  analysis  after  reduction.  2,848  were  type  I  vehicles,  and 
1,565  were  type  II.  Type  III  vehicles  were  not  included  in  the  analysis  because  their  CO 
and  HCs  exhaust  emission  concentrations  were  not  measured  due  to  the  fact  that  their 
sampling  points  were  in  a  different  plane.  The  analysis  of  the  retained  data  showed  that 
the  instantaneous  CO  and  HCs  exhaust  emission  concentrations  were  very  variable  and 
their  distributions  severely  skewed.  For  example,  the  average,  median,  standard 
deviation,  minimum,  and  maximum  CO  exhaust  emission  concentrations  for  the 
combined  fleet  were  1.12  %,  0.64%,  1.61  %,  0.01  %,  and  19.69  %  respectively.  Those  of 
HCs  exhaust  emission  concentration  for  the  combined  fleet  were  0.23  %,  0.20  %,  0.17  %, 
0.00%,  and  1.40  %  respectively.  The  mean  vehicle  speed  and  acceleration  were  56.56 
mph  and  0.17  mph/sec  respectively  for  the  combined  fleet.  Setting  high  emitters 
concentration  cut-points  to  4.0  %  CO  and  0.2  %  CO,  it  was  found  that  5.7%  and  49.2%  of 
the  combined  vehicle  fleet  were  CO  and  HCs  high  emitters  respectively. 

The  PVEMP  database  was  then  used  to  develop  a  modal  CO  and  HCs  exhaust 
emission  concentration  model.  The  underlying  theory  used  a  multivariate  regressive 
approach  applied  to  the  speed-acceleration  matrices  of  the  modal-average  CO  and  HCs 
exhaust  emission  concentrations  to  determine  the  three-dimensional  mathematical 
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relationships  between  mode-average  CO  and  HCs  exhaust  emission  concentrations  and 
speed  and  acceleration  midpoints.  The  regressions  were  found  to  fit  the  data  very  well. 
For  example,  the  multiple  coefficient  of  determination  was  equal  to  0.9792  for  the 
combined  fleet,  0.9934  for  type  I  vehicles,  and  0.9682  for  type  II  vehicles.  Based  on  the 
mathematical  relationships,  the  real-time  CO  and  HCs  exhaust  emission  concentration 
model  (PVEM-BEA)  was  developed  in  FORTRAN  for  the  PC  computing  environment. 
It  may  be  executed  as  an  integrated  program  to  the  Borman  Expressway  Traffic 
Management  Center  or  as  a  stand-alone  mobile  program.  It  requires  real-time  input 
parameters  (speeds,  vehicle  class,  and  road  geometry)  from  advanced  traffic  monitoring 
systems  such  as  the  Autoscope™.  It  can  also  get  its  inputs  from  two  consecutive  pre- 
installed  conventional  loops  as  long  as  means  to  assess  real-time  vehicle  types  exist. 
The  implementation  of  the  research  is  the  integration  of  the  real-time  modal 
model  to  the  Borman  Expressway  ATMS.  It  represents  a  tool  to  estimate  real-time 
vehicular  exhaust  emission  concentration,  especially  in  construction  work  zones. 
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APPENDIX  A 
A  Portion  of  the  Integrated  PVEMP  Database 
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Date 

rime 

CO  (%) 

HC  (%)         VT 

VS  (mph) 

VA(mph/s)T(oF) 

Td 

(OF) 

WS(mph)   WD(o)          RH(%) 

Sky 

5/19/97 

13:16:21 

1.2427110 

0.3569540 

1 

54.0 

-1.0 

57 

49 

10 

320 

74  OVC 

5/19/97 

13:16:22 

2.0345000 

0.3265460 

1 

54.0 

-1.0 

57 

49 

10 

320 

74  OVC 

5/19/97 

13:16:24 

7.2548070 

0.1647530 

1 

54.0 

-1.0 

57 

49 

10 

320 

74  OVC 

5/19/97 

13:16:26 

0.1753420 

0.2539020 

2 

54.0 

-1.0 

57 

49 

10 

320 

74  OVC 

5/19/97 

13:16:27 

0.5815650 

0.1694290 

2 

54.0 

-1.0 

57 

49 

10 

320 

74  OVC 

5/19/97 

13:16:29 

0.3967560 

0.3135310 

2 

52.0 

00 

57 

49 

10 

320 

74  OVC 

5/19/97 

13:16:31 

0.8495580 

0.1339370 

2 

52.0 

0.0 

57 

49 

10 

320 

74  OVC 

5/19/97 

13:16:32 

0.3557870 

0.2803900 

2 

52.0 

0.0 

57 

49 

10 

320 

74  OVC 

5/19/97 

13:16:34 
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Figure  A.l.  A  portion  of  the  integrated  PVEMP  database. 
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APPENDIX  B 
Graphical  Representations  for  Data  Pertaining  to  the  Two  Vehicle  Types 
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Figure  B  1.1.1.  Histogram  of  the  CO  exhaust  emission  concentration  distribution  for  type 
I  vehicles. 
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Figure  Bl.1.2.  Histogram  of  the  HCs  exhaust  emission  concentration  distribution  for  type 
I  vehicles. 
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Figure  B  1.2.1.  Histogram  of  the  CO  exhaust  emission  concentration  distribution  for  type 
II  vehicles. 
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Figure  Bl.2.2.  Histogram  of  the  HCs  exhaust  emission  concentration  distribution  for  type 
II  vehicles. 
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Figure  B2.1.1.  Logarithmic  distribution  of  CO  exhaust  emission  concentration  for  type  I 
vehicles. 
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Figure  B2.1.2.  Logarithmic  distribution  of  HCs  exhaust  emission  concentration  for  type  I 
vehicles. 
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Figure  B2.2.1.  Logarithmic  distribution  of  CO  exhaust  emission  concentration  for  type  II 
vehicles. 
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Figure  B2.2.2.  Logarithmic  distribution  of  HCs  exhaust  emission  concentration  for  type 
II  vehicles. 
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Figure  B3.1.  CO  and  HCs  exhaust  emission  concentration  percentiles  for  type  I  vehicles. 


1.6 


1.4 


.2    1.0 
is 

£    0.8 

0) 

o 

i    0.6 
o 

g    0.4 


0.2 


0.0 


CO  and  HC  Concentration  Percentiles 
Type  II  vehicles 


J. 

HC 

.■ 

" 

— y—  CO 

■• 

- 

■• 

.. 

" 

■■ 

■■ 

■■ 

■• 

■■ 

■■ 

•■ 

■■ 

■±-i^-rc~\~r 

i   i  i   i  l  i   i  i 

i  i   i   l  i   i 

Y  i  i  i  i  i 

i  i  i  i  i  i 

Mil- 

20 


15 


c 
o 


10  ■£ 
o 
u 
c 
o 
o 

5  ° 
0  O 


0 


10   20 


30 


40    50    60 
Percentiles 


70 


80 


90   100 


Figure  B3.2.  CO  and  HCs  exhaust  emission  concentration  percentiles  for  type  II  vehicles. 
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Figure  B4.1.1.  Variation  of  instantaneous  CO  exhaust  emission  concentrations  with 
vehicle  speed  for  type  I  vehicles. 
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Figure  B4.1.2.  Variation  of  instantaneous  CO  exhaust  emission  concentrations  with 
vehicle  acceleration  for  type  I  vehicles. 
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Figure  B4.1.3.  Variation  of  instantaneous  HCs  exhaust  emission  concentrations  with 
vehicle  speed  for  type  I  vehicles. 
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Figure  B4.1.4.  Variation  of  instantaneous  HCs  exhaust  emission  concentrations  with 
vehicle  acceleration  for  type  I  vehicles. 
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Figure  B4.2.1.  Variation  of  instantaneous  CO  exhaust  emission  concentrations  with 
vehicle  speed  for  type  II  vehicles. 
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Figure  B4.2.2.  Variation  of  instantaneous  CO  exhaust  emission  concentrations  with 
vehicle  acceleration  for  type  II  vehicles. 
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Figure  B4.2.3.  Variation  of  instantaneous  HCs  exhaust  emission  concentrations  with 
vehicle  speed  for  type  II  vehicles. 
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Figure  B4.2.4.  Variation  of  instantaneous  HCs  exhaust  emission  concentrations  with 
vehicle  acceleration  for  type  II  vehicles. 
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Figure  B5.1.1.  Three-dimensional  plot  for  the  triple  correlation  of  CO  emission 
concentration,  speed,  and  acceleration  for  type  I  vehicles. 
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Figure  B5.1.2.  Three-dimensional  plot  for  the  triple  correlation  of  HCs  emission 
concentration,  speed,  and  acceleration  for  type  I  vehicles. 
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Figure  B5.2.1.  Three-dimensional  plot  for  the  triple  correlation  of  CO  emission 
concentration,  speed,  and  acceleration  for  type  II  vehicles. 
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Figure  B5.2.2.  Three-dimensional  plot  for  the  triple  correlation  of  HCs  emission 
concentration,  speed,  and  acceleration  for  type  II  vehicles. 
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APPENDIX  C 
Graphical  Representations  of  the  Regression  Data  for  the  Two  Vehicle  Types 
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Figure  CI.  1.1.  Bin-average  CO  exhaust  emission  concentration  versus  speed  midpoints 
for  type  I  vehicles. 
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Figure  CI.  1.2.  Bin-average  HCs  exhaust  emission  concentration  versus  speed  midpoints 
for  type  I  vehicles. 


77 


o 
o 


u 

9   - 

A 

vg  CO    Type  1 

8   - 

7    - 

6   - 

5   - 

4    - 

• 

3   - 

2    - 

• 

• 
• 

1    - 
n 

• 

u      i 

1       '       i       •       i       ■       i 

1       ,       |       ,       |       ,       |       ,       | 

-3 


-2-10  1  2 

Acceleration  (mph/s) 


Figure  CI. 2.1.  Bin-average  CO  exhaust  emission  concentration  versus  acceleration 
midpoints  for  type  I  vehicles. 
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Figure  CI. 2. 2.  Bin-average  HCs  exhaust  emission  concentration  versus  acceleration 
midpoints  for  type  I  vehicles. 
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Figure  C2.1.1.  Bin-average  CO  exhaust  emission  concentration  versus  speed  midpoints 
for  type  II  vehicles. 
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Figure  C2.1.2.  Bin-average  HCs  exhaust  emission  concentration  versus  speed  midpoints 
for  type  II  vehicles. 
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Figure  C2.2.1.  Bin-average  CO  exhaust  emission  concentration  versus  acceleration 
midpoints  for  type  II  vehicles. 
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Figure  C2.2.2.  Bin-average  HCs  exhaust  emission  concentration  versus  acceleration 
midpoints  for  type  II  vehicles. 
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3D  CO  Type  I  mode-avg 


Figure  C3.1.1.  Three-dimensional  plot  of  the  mode-average  exhaust  CO  emission 
concentration  versus  speed  and  acceleration  midpoints  for  type  I  vehicles. 
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3D  HC  Type  I  mode-avg 


Figure  C3.1.2.  Three-dimensional  plot  of  the  mode-average  exhaust  HCs  emission 
concentration  versus  speed  and  acceleration  midpoints  for  type  I  vehicles. 
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3D  CO  Type  2  mode-avg 


Figure  C3.2.1.  Three-dimensional  plot  of  the  mode-average  exhaust  CO  emission 
concentration  versus  speed  and  acceleration  midpoints  for  type  II  vehicles. 
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3D  HC  Type  II  mode-avg 


Figure  C3.2.2.  Three-dimensional  plot  of  the  mode-average  exhaust  HCs  emission 
concentration  versus  speed  and  acceleration  midpoints  for  type  II  vehicles. 
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Figure  C4.1.1.  Three-dimensional  plot  of  the  mode-average  CO  exhaust  emission 
concentration  obtained  from  the  mathematical  relationship  for  type  I  vehicles. 
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Figure  C4.1.2.  Three-dimensional  plot  of  the  mode-average  HCs  exhaust  emission 
concentration  obtained  from  the  mathematical  relationship  for  type  I  vehicles. 
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Figure  C4.2.1.  Three-dimensional  plot  of  the  mode-average  CO  exhaust  emission 
concentration  obtained  from  the  mathematical  relationship  for  type  II  vehicles. 
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Figure  C4.2.2.  Three-dimensional  plot  of  the  mode-average  HCs  exhaust  emission 
concentration  obtained  from  the  mathematical  relationship  for  type  II  vehicles. 
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APPENDIX  D 
PVEM-BEA  Source  Code 
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program  pveml 
******************************************************************************** 

*The  Purdue  Micro-scale  Real-time  Vehicle  Exhaust  Emission  Concentration  Modal  * 
♦Model  version  1.0  (PVEM1.0) :  The  Borman  Expressway  Application  (PVEM-BEA  1.0) .* 
♦Developed  by  Dr.  Ouattara  Fatogoma,  Civil  Engineering,  Purdue  University,  * 
******************************************************************************** 

c 

C////MAIN  PROGRAM/////////////////////////////////////////////////////////////// 
c 

c Subroutines  definitions 

c  writerl,2:  output  description  of  the  model  and  all  the  modeling  data. 

c  readerl,2:  read  all  the  input  variables  and  calculates  the  #  of  vehicles. 

c  calcaccel :  calculates  average  speed  and  acceleration  for  each  vehicle. 

c  calcemission:  calculates  average  CO  and  HC  emission  concentations  for  each  vehicle 

c  and  each  vehicle  type. 

c  decmaker :  makes  decisions  about  emission  concentration  levels  per  vehicles  type. 

c  reseter:  resets  counters  and  variables  for  the  next  prediction  time. 

c  opener:  opens  all  the  necessary  files. 

c  closer:  closes  all  the  opened  files. 

c  initializer:  initializes  and  sets  all  constants  and  counters. 

c 

c Variables  definitions 

c  sp  =  average  vehicle  speed (mph)  within  averaging  time. 

c  msp  =  measured  vehicle  speed (mph) . 

c  ac  =  average  vehicle  acceleration (mph/s)  within  averaging  time. 

c  vt  =  vehicle  type ( typel (1) , type I I (2) , typelll (3) ) . 

c  em  =  average  exhaust  CO  and  HC  emission  concentrations (%) from  each  vehicle. 

c  aem  =  average  CO,HC  emission  concentrations (%)  per  vehicle  type. 

c  asp,aac  =  average  speed (mph)  and  acceleration (mph/s)  by  vehicle  type. 

c  eth  =  CO  and  HC  emission  concentration  thresholds ( %) . 

c  fig  =  CO  and  HC  emission  concentration  threshold  flags. 

c  eofg  =  end  of  file  flag. 

c  nv  =  number  of  vehicle  by  type . 

c  frq  =  frequency  of  vehicles  by  type  (#veh/s) . 

c  nl  =  number  of  traffic  lanes 

c  dl  =  distance  between  speed-reading  points  or  between  loops (ft) 

c  nvat  =  number  of  vehicles  emitting  more  than  the  thresholds. 

c  ts  =  time(hhmmss)  of  simulation. 

c  sss,mms,hhs  =  second,  minute,  hour  of  simulation. 

c  ds  =  date(mmddyy)  of  simulation. 

c  ti  =  time  interval  between  two  consecutive  speed  readings. 

c  eat  =  emission  concentration  averaging  time  (sec) . 

c  loc  =  name  of  the  emission  prediction  location. 

c  run  =  modeling  run  number. 

c  pef  =  product  of  emission  concentration  and  frequency  by  vehicle  type(%/sec) 

c  peft  =  product  of  emission  concentration  and  frequency  threshold (%/sec)  by  vehicle  typ 

e 

c  peffg  =  product  of  emission  concentration  and  frequency  flag 

c  dl  =  array  dimension  representing  number  of  vehicles  accross  speed  detection  points  in 

eat 
c  d2  =  array  dimension  representing  number  of  loops  in  eat 
c  x  =  array  dimension  representing  number  of  pollutants 
c  y  =  array  dimension  representing  of  vehicle  types 
c  ddm,tdm  =  date  and  time  of  decision  making 

c  cO, 1,2, 3,4, 5  =  mode-average  emission  concentration  regression  parameters 
c  bO, 1,2  =  accel  bin-average  emission  concentration  regression  parameters 
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c 

c Variables  declarations 

integer  dl,d2,x,y 

parameter (dl=100 , d2=2 , x=2 , y=4) 

integer  ds (dl, d2) , ts (dl,d2) , nl,msp (dl, d2) , nv(y) , run, nvat (x,y) , 
&vt (dl,d2) , rvt (dl) ,ddm, tdm, tnv 

real  eat, em(x,dl) , aem(x,y) , eth(x,y) , ti (dl) , asp (y) , aac (y) , dl, 
&sp(dl) , ae(dl) ,frq(y) ,pef (x,y) ,peft (x,y) , cO (x,y) , cl (x,y) ,c2 (x,y) , 
&c3 (x,y) ,c4 (x,y) , c5 (x,y) ,b0(x,y) ,bl(x,y) ,b2 (x,y) 

character  f lg(x,dl) *2, loc*5  0, eof g*4 , pef f g (x, y) *2, name*2  5 
c 

c Include  constants  initialization  file 

c      print  *,' Simulation  initialization  filename:  ' 

c      read (*, 10  0)  name 

c   100  format (a25) 

c      open (unit=9 , f ile=name) 

c       include  name 

include  ' c : \workpu\borman\programs\pvembeac . for ' 
c 
c Open  files 

call  opener 
c 
c Input  location  name, number  of  lanes, and  distance  between  loops 

call  readerl (loc, nl, dl) 
c 
c Output  model  description 

call  writerl (loc, nl, dl , eth, pef t) 
c 
c Input  vehicle  and  traffic  parameters.  Calculate  numbers  of  vehicles  and  frequencies 

10  call  reader2 (ds, ts, msp, vt , rvt , ti , eat , nv, nl , tnv, f rq, eof g, ddm, tdm) 
c 
c Calculate  average  speed  and  acceleration  for  each  vehicle  and  each  vehicle  type 

call  calcaccel (tnv, msp, sp, ac, ti , rvt , asp, aac, nv) 
c 

c Calculate  average  CO  and  HC  emission  concentrations  for  each  vehicle  and  and  each 

c    vehicle  type.  Flags  high  emitters.  Calculate  products  frequency  and  emission 

call  calcemission (tnv, nv, vt , em, bO , bl, b2 , cO , cl , c2 , c3 , c4 , c5 , eth, fig, 
&  nvat , sp, ac, aem, pef , frq,  nl) 

c 

c    Make  decission  about  vehicle  type  whose  emission  concentration  exceed  designed  thre 
sholds 

call  decmaker (pef , pef t , pef f g, aem, ddm, tdm) 
c 
c Output  individual  vehicle  modeling  data  and  total  fleet  modeling  data 

call  writer2 (tnv, run, ds, ts, rvt, sp, ac, em, f lg, nv, frq, nvat , 
&  asp, aac, aem, pef , pef fg) 

c 
c Reset  counters  for  the  next  simulation  time  or  end  the  simulation  if  end  of  file 

call  reseter (run, eofg, *10) 
c 
c End  main  program 

stop 

end 
c 
c 
C////SUBROUTINES  FOR  FILE  I/O//////////////////////////////////////////////////// 


91 

c Files  description 

c  pvembeai.dat  =  traffic  and  vehicle  parameters  input  file 

c  pvembeao.dat  =  modeling  output  file 

c  pvembeac.dat  =  simulation  constant  initialization  input  file 

c 

c Open  all  necessary  files 

subroutine  opener 
c 

c Simulation  input  and  output  files 

c      character  namei*24 , nameo*24 

c      print  *, 'Enter  traffic  data  input  filename:  ' 

c       read(*,10  0)  namei 

c      print  *, 'Enter  modeling  output  filename:  ' 

c      read (*, 10  0)  nameo 

c   100  format (a24) 

c      open (unit=10 , f ile=namei) 

c      open (unit=20 , f ile=nameo) 

open (unit=10 , f ile= ' c : \workpu\borman\programs\pvembeai . dat ' ) 

open (unit=20 , f ile= ' c : \workpu\borman\programs\pvembeao .dat ' ) 
c 

c Simulation  constants  file 

c 

c Position  all  files  to  the  beginning 

rewind  10 

rewind  2  0 
c 

end 
c 
c 

c Close  input  and  output  files 

c 

subroutine  closer 

close (10) 

close (20) 

end 
c 
c 

C////SUBROUTINE  FOR  VARIABLES  RESET///////////////////////////////////////////// 
c 

subroutine  reseter (rr, reof g, *) 
c 

c Variables  definition 

c  rr  =  Run  number  counter 

c  reofg  =  end  of  file  flag 

c 

c variable  declaration 

integer  rr 

character  ryon,reofg*4 
c 
c Next  modeling  session  w/o  interuption 

if (reof g. eq. ' fend ' )  then 
call  closer 

else 

rr=rr+l 
return  1 

endif 
c 
c Next  modeling  session  with  pause 
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c    30  print  *, 'Do  you  wish  to  run  another  simulation?  (y  or  n) :  ' 

c       read ( * , 1 0  0 )  ryon 

c       if (ryon.eq. 'y ' )  then 

c        if (reof g.ne . ' f end ' )  return  1 

c        call  closer 

c       elseif (ryon. eq. ' n' )  then 

c         call  closer 

c       else 

c        print  *,  'Enter  y  or  n!  ' 

c         goto  3  0 

c       endif 

c 

100  format (al) 

end 
c 
c 

C////SUBROUTINE  FOR  READING  INPUT  DATA////////////////////////////////////////// 
c 

c Read  location  name,  number  of  traffic  lanes,  distance  between  speed  detection  point 

s 
c 

subroutine  readerl (rloc, rnl , rdl) 
c 

c Variables  definition 

c  rloc  =  modeling  location 

c  rnl  =  number  of  traffic  lanes 

c  rdl  =  distance  between  the  2  loops 

c 

c Variables  declaration 

integer  rnl 

real  rdl 

character  rloc*50 
c 

c List  directed  input 

c      print  *, 'Location  name:  ' 

c       read (*, 110)  rloc 

c      print  *, 'Number  of  traffic  lanes:  ' 

c       read (*, 115) ,  rnl 

c      print  *, 'Distance  between  the  speed  detection  points:  ' 

c      read (*, 120) ,  rdl 

c 

c File  input 

read (10, 100)  rloc, rnl, rdl 
c 

100  format (a50/il/f 6 .2/) 
110  format (a50) 
115  format (il) 
120  format(f6.2) 

end 
c 
c 

c Read  all  vehicle  parameters:  vehicle  type,  speeds  at  loops 

c     Calculate  number  and  frequency  of  vehicles  by  type 
c 

subroutine  reader2 (rds, rts, rmsp, rvt , rrvt , rti, reat , rnv, rnl, rtnv, 
&  rf rq,  reof g, rddm, rtdm) 

c 
c Variables  definition 
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c  i , j  =  counters 

c  rts  =  time (hhmmss)  of  simulation. 

c  rsss, rmms, rhhs  =  second,  minute,  hour  of  simulation. 

c  rds  =  date (mmddyy)  of  simulation. 

c  rd  =  array  dimension  =  maximum  number  of  vehicles 

c  rsp  =  vehicle  speed (mph) 

c  rvt  =  vehicle  type 

c  rnv  =  car, truck, car+truck, heavy  trucks, combined  fleet  counts  at  loops 

c  ravt  =  averaging  time (sec) 

c  rtsec  =  time  in  second 

c  rti  =  time  interval  for  two  consecutive  readings ( sec) 

c  reat  =  maximum  averaging  time  interval (4.0  sec) 

c  reofg  =  end  of  file  flag 

c  rfrq  =  vehicle  frequency  by  type (#veh/sec) 

c  rddm,rtdm  =  date, time  of  decision  making 

c 

c Variables  declaration 

integer  rdl , rd2 , ry 

parameter  (rdl=100 , rd2=2 , ry=4) 

integer  i,  j  ,  1,  il,  i2 , rn, rmsp (rdl, rd2)  , rnv(ry) , rds (rdl, rd2)  , 
&rts (rdl,rd2) , rsss (rdl,rd2) , rmms (rdl,rd2) , rhhs (rdl, rd2) , 
&rnl, rvt (rdl , rd2) , rrvt (rdl) , rddm, rtdm, rtnv 

real  reat, rti (rdl) , rtsec (rdl, rd2) , rfrq (ry) 

character  ryon,reofg*4 
c 
c counter  initialization  and  constants 

i  =  0 

j-0 

rtnv=0 

do  5  m=l,4 
rnv(m) =0 
5  continue 


ccList  directed  input  of  real-time  traffic  data 

c First  loop  data  reading 

20  j=j+l 
30  i=i+l 

print  ' (a\) ■ , 'Date  of  vehicle 

read  * , rds (i , j ) 

print  ' (a\) ' , 'Time  of  vehicle 

read  * , rts (i , j ) 

print  ' (a\) ' , 'Vehicle, 

read  *,rmsp(i,j) 

print  *, 'Vehicle, 

print  *,'l=type  I 

read  * , rvt ( i , j ) 


i  , 


at  point 
at  point 
speed  at  point  ' , 


,  i ,  '  type  at  point  ' , j . 
2=type  II,  3=type  III' 


continue  reading 

18  print  ' (a\) ' , 'Another  Vehicle?  (y  or  n) 
read  * , ryon 
if (ryon.eq. 'y' )  then 

goto  3  0 
elseif (ryon.eq. ' n' )  then 

print  ' (a\) ' , 'Reading  over  at  loop  ' , 
else 

print  ' (a\) ' , 'Enter  y  or  n' 

goto  18 
endif 


1  (mmddyy) : 
'  (hhmmss) : 
(mph) :  ' 
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c 

c Second  loop  data  reading 

c       if(j.lt.2)  then 

c        i  =  0 

c        goto  2  0 

c      endif 

c 

c 

ccFile  input  of  real-time  traffic  data 

c 

40  j=j+l 
50  i=i+l 
c 

c Read  data  recorded  at  point  j 

60  read(10,110)  rds (i , j ) , rts (i , j ) , rmsp (i , j ) , rvt ( i , j ) 
c 
c Convert  time  in  second 

rsss (i, j ) =mod (rts (i, j ) ,100) 

rmms(i, j ) =mod ( (int (rts(i, j) /100) ) ,100) 

rhhs(i, j)=int (rts(i, j) /10  000) 

rtsec (i, j ) =rhhs (i, j ) *3600 . 0  +  rmms (i, j ) *6  0 . 0  +  rsss (i,  j  )  *1 . 0 
c 
c Next  vehicle  at  that  loop 

if(i.lt.rnl)  goto  50 
c      if ( (rtsec (i, j ) -rtsec (1, j )). le . reat)  then 
c        goto  5  0 
c      else 

c        rds (1,2) =rds (i, j) 
c        rts (1,2) =rts (i, j ) 
c        rmsp (1, 2) =rmsp (i, j ) 
c        rvt (l,2)=rvt (i, j) 
c        j=2 
c        i  =  l 
c        goto  5  0 
c      endif 
c 
c Total  number  of  vehicles  at  point  j 

if (j .eq.l)  il=i 

if ( j .eq.2)  i2=i 
c 
c Read  data  recorded  at  point  j+1 

if(j.lt.2)  then 
i  =  0 
goto  4  0 

endif 
c 

c Check  for  data  consistency  in  vehicle  type, number  and  licence  plate  number, 

c     and  time  and  calculate  time  interval  and  real  vehicle  counts 
c Case  I&III:  no  change  in  lanes 

if(i2. eq.il)  rn=il 
c Case  II:  convergence;  eliminate  far  end  vehicle 

if(i2. lt.il)  rn=i2 
c Case  IV:  divergence;  speed  1  of  the  added  veh=speed  of  the  closest  veh 

if (i2.gt.il)  rn=i2 
c 

do  70  i=l,rn 
c        if (rlpn(i, 2) .eq.rlpn(i, 1) )  then 
if (rvt (i,2) .eq.rvt (i,l) )  then 
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rrvt ( i ) =rvt ( i , 1 ) 

rti (i) =rtsec ( i, 2) -rtsec (i , 1) 

c Calculate  number  of  vehicles  by  type  at  loop  j.  1  =  1 (typel)  , 1=2 (type2) , 1  =  3 (type 

3) 

c         1=4 (combined  1+4  fleet) 
c  if  (rti (i) . le .reat)  then 

if (rrvt (i) .eq.l)  1=1 
if (rrvt (i) .eq.2)  1=2 
if (rrvt (i) .eq.3)  1=3 
rnv ( 1 ) =rnv ( 1 ) + 1 

if ( (1 .eq.l) .or. (1 .eq.2) )  rnv (4) =rnv(4) +1 
rtnv=rtnv+l 
c  endif 

endif 
70  continue 

c 

c End  of  file  clause 

if(eof(10))  then 
reof g= ' fend' 

endif 
c 
c Vehicle  flow  rate  (#  vehicle  within  the  maximum  averaging  time  per  lane) 

do  80  i=l,4 

rf rq(i) = (rnv(i) /reat) /real (rnl) 
8  0  continue 
c 

rddm=rds (1,1) 

rtdm=rts(l, 1) 
c 

110  format (2 (is, lx) ,i2,lx,il) 
99  end 
c 
c 

C/// /SUBROUTINE  FOR  THE  CALCULATION  OF  VEHICLE  AVERAGE  SPEED  AND  ACCELERATION// 
c  calculates  the  average  speed  and  acceleration  for  each  vehicle  and  each 
c  vehicle  type  within  the  time  interval 
c 

subroutine  calcaccel (ctnv, cmsp, csp, cac, cti , cvt , casp, caac, cnv) 
c 

c Variables  definitions 

c  cd  =  arrays  dimension 

c  cmsp  =  measured  vehicle  speed (mph) 

c  cvt  =  vehicle  type 

c  csp  =  average  vehicle  speed  within  time  interval (mph) 

c  cac  =  average  vehicle  acceleration  within  time  interval (mph/sec) 

c  cssp  =  sum  of  speed  by  vehicle  type (mph) 

c  casp  =  average  speed  by  vehicle  type (mph) 

c  csac  =  sum  of  acceleration  by  vehicle  type (mph/sec) 

c  caap  =  average  acceleration  by  vehicle  type (mph) 

c  cti  =  time  interval (sec) 

c  ctnv  =  total  number  of  vehicles 

c  cnv  =  number  of  vehicles  by  type 

c 

c Variable  declaration 

integer  cdl,cd2,cy 

parameter (cdl=100 , cd2=2 , cy=4 ) 

integer  cmsp ( cdl , cd2 ) , i , j , 1 , k , cnv ( cy ) , cvt ( cdl ) , ctnv 
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c 
c- 


real  csp (cdl) , cac (cdl) , cti (cdl) , cssp (cy) , csac (cy) , casp (cy) , 
&caac (cy) 

—  Initialization 

do  5  1=1,4 
cssp(l)  =0  .  0 
csac(l) =0 . 0 
5  continue 

—  Calculate  average  speed  and  acceleration  within  time  interval 

do  10  i=l,ctnv 

csp  (i)  =real  (cmsp  (i ,  2)  +cmsp  (i,  1)  )  /2  .  0 
cac (i)  =real (cmsp (i, 2) -cmsp (i,l))/cti(i) 


-Range  accel 
if (cac (i) 


if  ( 

if  ( 
if  ( 
if  ( 
if  ( 
if  ( 
if  ( 
if  ( 
if  ( 
if  ( 
if  ( 
if  ( 
if  ( 
if  ( 
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cac (i) =-1 , 

cac  (i) =-0 . 
cac (i) =0 . 0 
cac (i) =0  .  5 
cac (i) =1 . 
cac ( i) =1 . 
cac (i) =2 . 
cac (i) =2 . 
cac  (  i) =3 . 
cac ( i) =3 . 


c 
c- 


c 

c- 


Sum  of  speeds  and  accelerations  by  vehicle  type 

if (cvt (i) .eq.l)  j=l 

if (cvt (i) .eq.2)  j=2 

if (cvt (i) .eq.3)  j=3 

cssp ( j ) =cssp ( j ) +csp  (i) 

csac ( j ) =csac ( j ) +cac (i) 

if ( ( j .eq.l) .or . ( j .eq.2) )  then 

cssp(4)=cssp(4)+csp(i) 

csac (4) =csac (4 ) +cac (i) 
endif 
10  continue 

— Calculate  average  speed  and  acceleration  by  vehicle  type 
do  15  k=l,4 

if (cnv(k) .eq. 0)  then 
casp(k) =0 . 0 
caac (k) =0 . 0 
else 

casp (k) =cssp (k) /real (cnv(k)  ) 
caac (k) =csac (k) /real (cnv(k) ) 
endif 
15  continue 
end 
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C////SUBROUTINE  FOR  THE  CALCULATION  OF  VEHICLE  MODE-AVERAGE  EMISSION  CONCENTRATIONS/ 

c  Calculates  the  average  co  and  he  emission  concentrations  for  each  vehicle  and 

c  eack  vehicle  type  as  well  as  the  product  of  concentration  and  frequecy.  Flags 

c  all  the  high  emitters 

c 

subroutine  calcemission (ctnv, env, cvt ,  cem, cbO, cbl , cb2 ,cc0,ccl,cc2, 
&  cc3,cc4,cc5,ceth,cflg, envat , csp, cac, caem, 

&  cpef , cf rq, cnl) 

c 

c Variables  definitions 

c  cd,cx,cy  =  arrays  dimensions 

c  i,j,k,l,m,n  =  counters 

c  en  =  number  of  vehicles  by  type  for  emission  averaging  purpose 

c  envat  =  number  of  vehicles  by  type  above  thresholds 

c  csp  =  vehicle  speed (mph) 

c  cac  =  vehicle  acceleration (mph/s) 

c  cvt  =  vehicle  type 

c  cem  =  CO  and  HC  average  emission  concentrations  within  time  interval (%) 

c  ceth  =  CO  and  HC  emission  concentration  thresholds (%) 

c  csem  =  sum  of  emission  concentrations  by  vehicle  type  w/i  time  interval (%) 

c  caem  =  average  emission  concentrations  by  vehicle  type  w/i  time  interval (%) 

c  cflg  =  emission  concentration  threshold  flags 

c  cpef  =  product  of  emission  concentration  and  vehicle  frequency  by  type(%/s) 

c  cfrq  =  vehicle  frequency (#veh/s) 

c  ccO, 1, 2 , 3 , 4 , 5  =  mode-average  emission  concentration  regression  parameters 

c  cb0,l,2  =  accel  bin-average  emission  concentration  regression  parameters 

c 

c Variables  declaration 

integer  cdl , ex, cy, i , k, 1 , m, n 
parameter (cdl=l 00, cx=2 , cy=4) 

integer  en (ex, cy) , envat (ex, cy) , cvt (cdl) , env (cy) , cnl , ctnv 
real  csp (cdl) , cac (cdl) , cem (ex, cdl) , ceth (ex, cy) , csem (ex, cy) , 
&caem (ex, cy) , cpef (ex, cy) , cf rq (cy) , ccO (ex, cy) , ecl (ex, cy) , cc2 (ex, cy) , 
&cc3  (ex,  cy)  ,  cc4  (ex,  cy)  ,  cc5  (ex,  cy)  ,  cbO  (ex,  cy)  ,  cbl  (ex,  cy)  ,  cb2  (ex,  cy) 
character  cflg (ex, cdl) *2 
c 

c Constants  initializations 

k=0 
c 

c Emission  summation, number, high  emitter  flags  initialization 

do  10  m=l,2 
do  15  n=l,4 

csem (m, n) =0 . 0 
en (m, n) =0 
envat (m, n) =0 
15    continue 
10  continue 
c 

c Calculate  emission  concentrations 

c    k=l  for  CO  emission  concentration,  k=2  for  HC  emission  concentration 
c     1=1  for  type  1,1=2  for  type  2,1=3  for  type  3,1=4  for  combined  fleet 
c 

20  k=k+l 

do  3  0  i=l,ctnv 

if (cvt (i) .eq.l)  1=1 
if (cvt (i) .eq.2)  1=2 
if (cvt(i) .eq.3)  1=3 
c Vehicle  numbers  for  emission  averaging 
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cn(k/l)=cn(k,l) +1 

if ( (l.eq.l) .or. (l.eq.2) )  en (k, 4 ) =cn (k, 4 ) +1 
c 

c Emission  calculation  formula 

if (csp(i)  .It .35.  0)  then 
cem(k, i) =cb0 (k, 1) + 
&       cbl(k,l)  *(  (cac(i) -0.0)/4.0)  + 
&        cb2 (k, 1) * ( (cac(i) -0 . 0) /4 . 0) **2 . 0 
else 

cem (k,i)=cc0(k,l)+ 
&        ccl(k,l) * ( (csp(i) -55.0)/30.0) + 
&        cc2 (k, 1) * ( (csp(i) -55.0)/3  0.0)**2.0+ 
&        cc3 (k, 1) * ( (cac(i) -0. 0) /4 . 0) + 
&        cc4 (k,l) * ( (cac(i) -0.0)/4.0)**2.0+ 

&        cc5(k,l) * ( (csp(i) -55.0)/3  0.0)*( (cac(i) -0.0) /4.0) 
endif 
c 

c High  emitters.  at=above  emission  concentration  threshold,  bt=below. 

if (cem (k, i) .gt . ceth (k, 1) )  then 
cflg(k,i)='at' 
envat (k,  1) =cnvat (k, 1) +1 
else 

cflg(k,i)='bt ' 
endif 
c 

c Sum  of  emission  concentrations 

csem(k, 1) =csem(k, 1) +cem (k, i) 

if ( (1 .eq. 1) .or. (1 .eq.2) )  csem (k,4)=csem(k,4)+cem(k,i) 
c 

3  0  continue 
c 

c Next  pollutant 

if(k.ne.2)  goto  20 
c 

c Average  emission  concentration  and  product  of  (average  emission/number  of  vehicle) 

and 
c      flow  rate  and  number  of  lane  per  vehicle  type 
do  4  0  m=l,2 
do  5  0  n=l,4 

if (cn(m,n) .eq.0)  then 

caem (m, n) =0 . 0 
else 

caem (m,n) =csem (m, n) /real (cn(m,n) ) 
endif 

cpef (m,n) =caem(m, n) *cf rq (n) *cnl 
50    continue 

4  0  continue 

end 
c 
c 

C////SUBROUTINE  FOR  DECISION  MAKING///////////////////////////////////////////// 
c  Flags  the  vehicle  types  which  average  emission  concentrations  by  frequency 
c  products  are  greater  than  the  designed  Thresholds 
c 

subroutine  decmaker (dpef , dpef t , dpef f g, daem, dtdm, dddm) 
c 

c Varibles  definitions 

c  dpef  =  product  of  emission  concentration  and  vehicle  frequency (%/sec) 
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c  dpeft  =product  of  emission  concentration  and  vehicle  frequency  threshold ( %/sec) 
c  dpeffg  =  product  of  emission  concentration  and  vehicle  frequency  flag 
c  daem  =  average  emission  concentration  by  vehicle  type 
c 

c Variables  declaration 

integer  i, j , dx, dy, dtdm, dddm 
parameter (dx=2 , dy=4 ) 

real  dpef (dx, dy) , dpeft (dx, dy) , daem (dx, dy) 
character  dpef f g (dx, dy) *2 
c 

do  10  i-1,2 
do  20  j=l,4 
dpeffg (i, j)='bt' 

if (dpef (i, j) .gt .dpeft (i, j) )  then 
dpef fg(i, j ) = ' at ■ 

write ( * , *)  dtdm, dddm, i , j , daem (i , j ) 
call  beepqq( 1000, 150) 
endif 
20    continue 
10  continue 
end 


C////SUBROUTINES  FOR  MODELING  RESULTS  AND  MODEL  DESCRIPTION  OUTPUT////////////// 

c 

c 

c Ouput  aknowledgements  and  descriptions,  and  simulation  constants 

c 

subroutine  writerl (wloc, wnl, wdl , weth, wpef t ) 
c 

c Variables  definition 

c  wloc  =  location  name 

c  wnl  =  number  of  traffic  lanes 

c  wdl  =  distance  between  traffic  lanes (ft) 

c  weth  =  indivudual  vehicle  emission  concentration  threshold (%) 

c  wpeft  =  product  of  emission  concentration  and  vehicle  frequency  threshold (%/sec) 

c 

c Variable  declaration 

integer  wx,wy,j,wnl 

parameter  (wx=2,wy=4) 

real  wdl , weth (wx, wy) , wpeft (wx, wy) 

character  wloc*50 
c 

write(20, 95) 

write(20,210) 

write(20, 97) 

write(20,101) 

do  5  j=l,4 

write (2  0, 103)  j , weth(l, j ) , weth (2, j ) 
5   continue 

write(20, 102) 

do  10  j=l,4 

write (2  0, 103)  j , wpeft (1, j ) , wpef t (2, j ) 
10   continue 

write (20, 200) 

write (20, 100)  wloc, wnl, wdl 

write(20,210) 


100 


95  format ( '  ******************************************************* • 
£  f  ********************* i  / 

&'  *  Purdue  Real-Time  Vehicle  Exhaust  Emission  Modal  Model  V  1.0', 
&lx, ' (PVEM1 . 0) ' ,4x, ' * ' / '  *  Developed  by  Dr.  Ouattara  Fatogoma.  ', 
&' Civil  Engineering,  Purdue  University  *'/ 
&'  *  Copyright  (C)  1998.  All  Rights  Reserved .', 33x, '*' / 

r,  I   **************************************************************! 
&  I  **************  l  //) 

97  format ( 'NOTES '/'=====' ) 

100  format (/ 'Modeling  location:  'a50  /'Number  of  lanes:  'il/ 
^'Distance  between  2  speed  detection  points  (ft) :  'f6.2) 

101  f ormat ( ' CO  and  HCs  instantaneous  exhaust'/ 
^'emission  concentration  thresholds (%) '/ 
&'VT' ,4x, 'CO' , 6x, 'HCs' ) 

102  format (/ 'Designed  CO  and  HCs  threshold  limit  values'/ 
&'by  vehicle  type(%/sec) '/ 

Sc'VT'  ,4x,  'CO'  ,6x,  'HCs'  ) 

103  format (il,2x, f 7 .4,2x, f7 .4) 

200  format (/ 'Notes :  bt  =  below  threshold,  at  =  above  threshold'/ 
fc'l=type  1  vehicles  (LDV) ,  2=type  2  vehicles  (MDV) , ' , lx, 
&'3=type  3  vehicles  (HDV) ' ,  '4=Combined  1&2 ' ) 

210  format ( '=======================================================' , 

&.=====================================./) 

end 
c 
c 
c Output  vehicle,  traffic,  and  emission  data  subroutine 

subroutine  writer2 (wtnv, wrun, wds, wts, wvt , wsp, wac , wem, wf lg, wnv, 
&  wf rq, wnvat , wasp, waac, waem, wpef , wpef f g) 

c 

c Variables  definition 

c  wrun  =  simulation  run  number 

c  wds  =   date  of  simulation (mmddyy) 

c  wts  =   time  of  simulation (hhmmss) 

c  wvt  =   vehicle  type 

c  wsp  =   vehicle  speed (mph) 

c  wac  =   vehicle  acceleration (mph/s) 

c  wem  =   average  emission  concentration (%) 

c  wflg  =  emission  threshhold  flag 

c  wnv  =   vehicle  count  by  type 

c  wfrq  =  vehicle  frequency (%/s) 

c  wnvat  =  number  of  vehicles  above  threshold 

c  wasp  =   average  vehicle  speed  by  type (mph) 

c  waac  =   average  vehicle  acceleration  by  type (mph/s) 

c  waem  =   average  emission  concentration  by  type(%) 

c  wpef  = 

c  wpeffg  = 

c 

c Variable  declaration 

integer  wdl, wd2 , wx, wy, i 

parameter  (wdl=100 , wd2=2 , wx=2 , wy=4) 

integer  wts(wdl,wd2) , wds (wdl , wd2) , wrun, wnv (wy) , wnvat (wx, wy) , 
Siwvt  (wdl)  ,wtnv 

real  wsp (wdl) , wac (wdl) , wem (wx, wdl) , wasp (wy) , waac (wy) , 
&waem (wx, wy) ,wfrq(wy) , wpef (wx,wy) 

character  wflg (wx, wdl) *2 , wpeffg (wx, wy) *2 
c 
c Write 
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,  wac  (i)  ,  wetn  (1 ,  i)  , 


output  per  type 
wnv ( 1 ) , wnv ( 2 ) , wnv ( 3 ) , wnv ( 4 ) 
wfrq(l) ,wfrq(2) ,wfrq(3) ,wfrq(4) 
wasp(l) ,wasp(2) ,wasp(3) ,wasp(4) 
waac (1) , waac (2) , waac (3) , waac (4) 
waem (1,1) , waem (1,2) , waem (1,3) , waem ( 1 , 4 ) 
waem (2,1) , waem (2,2) , waem (2,3) , waem (2,4) 
wnvat (1,1) , wnvat (1,2) , wnvat (1,3) , wnvat (1,4) 
wnvat (2,1) , wnvat (2,2) , wnvat (2,3) , wnvat (2,4) 
wpef (1,1) ,wpef (1,2) ,wpef (1,3) ,wpef (1,4) 
wpef (2,1)  ,wpef (2,2)  , wpef (2,3)  , wpef (2,4) 
wpef fg(l, 1)  , wpef fg(l, 2)  , wpef fg(l, 3)  , wpef fg(l, 4) 
wpef fg(2, 1) , wpef fg(2, 2) , wpef fg(2, 3) , wpef fg(2, 4) 


write (20, 100)  wrun 

write(20,105) 

do  10  i=l,wtnv 

Write (2  0, 110)  wds (i,l) ,wts(i,l) , wvt (i) ,wsp(i 
&  wf lg(l, i) , wem(2, i) , wflg(2, i) 

10  continue 

c Output  header 

write(20, 120) 
write(20,121) 

c Average  data 

write(20, 130) 
write (20, 132) 
write (20, 135) 
write (20, 140) 
write(20,145) 
write (20, 150) 
write(20,155) 
write(20,160) 
write(20,165) 
write(20,170) 
write(20,175) 
write (20, 180) 
write (20,210) 

c Formats 

10  0  f ormat ( ' Run# : 

105  format ( 'SIMULATION  DATA '/'==============='/'  Date ', 5x, ' Time ' 

&'Vt ' ,3x, ' Spd' ,4x, 'Ace' , 7x, 'CO' ,4x, 'COfg' ,3x, ' HC ' ,4x, 

&'HCfg' /' (mmddyy) ' , ' (hhmmss) ' , 6x, ' (mph) ' , lx, ■ (mph/s) ' , 5x, '(%)', 

&10x, '(%)•) 

2 (i6 .6,3x) ,il,3x,2 (f5.2,3x) , f 7 . 4 , 2x, a2 , 2x, f 7 . 4 , 2x, a2) 
//'SUMMARY  DATA' /'=  =  =  =  =  =  =  =  =  =  =  =' ) 

44x, 'Type_l' ,4x, ' Type_2 ' ,4x, ' Type_3 ' ,4x, ' CF_1&2 ' ) 

'Number  of  vehicles:  ' , t44 , 4 (i3 , 7x) ) 

'Vehicle  Flow  rate (#veh/sec/lane) :  ' , t45 , 4 (f 5 . 2 , 5x) ) 

'Average  speed(mph):  '  , t45 , 4  (f 5 . 2 , 5x) ) 

'Average  acceleration (mph/sec) :  ' , t45 , 4 (f 5 . 2 , 5x) ) 

'Average  CO  emission  concentration (%) :  ' , t45 , 4 (f 7 . 4 , 3x) ) 

'Average  HC  emission  concentration (%) :  ' , t45 , 4 (f 7 . 4 , 3x) ) 

'Number  of  vehicles  exceeding  CO  threshold:  ' , t44 , 4 (i3 , 7x) ) 

'Number  of  vehicles  exceeding  HC  threshold:  ' , t44 , 4 (i3 , 7x) ) 

'CO  PCFNL(%/sec) :  '  , t45 , 4  (f 7 . 4 , 3x) ) 

'HCs  PCFNL(%/sec) :  ' , t45 , 4 (f 7 . 4 , 3x) ) 

'CO  PCFNL  flag:  '  , t46 , 4 (a2 , 8x) ) 

•HCs  PCFNL  flag:  ' , t46 , 4 (a2 , 8x) ) 

*  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =:  =  =  =  =  =  =  =  =  =:  =  =  =  ■ 


,i6.6/) 


,4x. 
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format ( 

120 

format ( 

121 

format ( 

130 

format ( 

132 

format ( 

135 

format ( 

140 

format ( 

145 

format ( 

150 

format ( 

155 

format ( 

160 

format ( 

165 

format ( 

170 

format ( 

175 

format ( 

180 

format  ( 

210 

format  ( 

end 
C///////////////////////////////////// 


